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The  Consolidated  Compartment  Fire  Model  (CCFM)  Computer  Code 
Application  CCFM. VENTS  - Part  I:  Physical  Basis 


Leonard  Y.  Cooper  and  Glenn  P.  Forney 


ABSTRACT 

A project  was  carried  out  at  The  National  Institute  of  Standards  and 
Technology  (NIST)  to  study  the  feasibility  of  developing  a new- generation , 
multi-room,  compartment  fire  model  computer  code,  called  the  Consolidated 
Compartment  Fire  Model  (CCFM)  computer  code.  The  idea  was  that  such  a code 
would  consolidate  past  progress  in  zone -type  compartment  fire  modeling,  and 
allow  readily  for  integration  of  future  advances  with  the  greatest  possible 
flexibility.  Desired  features  of  the  CCFM  would  include:  comprehensive 
documentation,  user-friendliness,  significant  modularity,  nrimerical 
robustness,  and  versatility  in  the  sense  that  the  code  would  provide  a 
capability  of  analyzing  a particular  compartment  fire  problem  by  using  any  one 
of  a range  of  physical-phenomena-modeling  sophistication,  from  the  most  basic 
to  the  most  comprehensive.  The  project  led  to  the  development  of  a prototype 
multi-room  CCFM  product  called  CCFM. VENTS.  CCFM. VENTS  involves  a model 
formulation  and  code  structure  that  allows  for  the  required  future  CCFM  growth 
flexibility.  It  has  a relatively  sophisticated  and  very  general  room-to-room 
forced  and  unforced  vent  flow  capability.  Finally,  the  CCFM. VENTS  code  uses 
the  simplest  possible,  point-source-plume,  smoke  - filling  fire  physics  in  the 
rooms -of - fire-origin  and  a very  simple  heat  transfer  calculation  there  and  in 
other  spaces. 

This  is  Part  I of  a four-part  report  which  documents  the  above  effort. 
Introductory  remarks  discuss  the  generic  features  of  the  CCFM  and  the  specific 
features  of  CCFM. VENTS.  The  main  objective  of  this  Part  I document  is  to 
present  a comprehensive  description  of  the  governing  equations  of  CCFM. VENTS 
and  their  technical  basis. 

The  other  three  parts  of  this  report  are:  Part  II:  Software  Reference  Guide; 
Part  III:  Catalog  of  Algorithms  and  Subroutines;  and  Part  IV:  User  Reference 
Guide . 

Some  of  the  features  of  CCFM. VENTS  are  now  in  the  process  of  being  integrated 
into  CFAST,  a revised  version  of  the  FAST  fire  model  which  is  being  used 
currently  in  the  NIST/CFR  Fire  Hazard  Assessment  Method,  HAZARD. 


Keywords : building  fires ; compartment  fires ; computer  models ; fire  models ; 

mathematical  models;  vents;  zone  models. 
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1 . Background 


1.1.  The  Consolidated  Compartment  Fire  Model  (CCFM) ; A New  Generation  Multi- 
Room  Fire  Model 

A project  was  carried  out  at  The  National  Institute  of  Standards  and 
Technology  (NIST)  to  study  the  feasibility  of  developing  a new- generation , 
multi-room,  compartment  fire  model  computer  code,  called  the  Consolidated 
Compartment  Fire  Model  (CCFM)  computer  code.  The  idea  was  that  such  a code 
would  consolidate  past  progress  in  zone- type  compartment  fire  modeling,  and 
allow  readily  for  integration  of  future  advances  with  the  greatest  possible 
flexibility.  The  CCFM  would  be  developed  in  stages  and  the  development 
process  would  lead  to  a series  of  well - documented , user-friendly,  numerically- 
robust,  modular,  and  easily  updated  CCFM  products  or  "applications,".  These 
would  be  versatile  in  the  sense  that  they  would  provide  a capability  of 
simulating  a particular  compartment  fire  scenario  by  using  one  of  a series  of 
available  levels  of  physical  modeling  sophistication  from  the  most  basic  to 
the  most  comprehensive.  The  final  stage  CCFM  would  provide  the  best  possible 
mathematical  simulation  of  fire  - generated  environments. 

The  project  led  to  a prototype  multi-room  CCFM  product  called  CCFM. VENTS 
(i.e.,  an  application  of  the  CCFM  concept  named  VENTS).  CCFM. VENTS  involves  a 
model  formulation  and  code  structure  that  allows  for  the  required  future  CCFM 
growth  flexibility. 

This  document.  Part  I;  Physical  Basis,  is  the  first  of  a four-part  report 
which  documents  CCFM. VENTS.  The  other  three  parts  are;  Part  II;  Software 
Reference  Guide  [1.1.1]^;  Part  III;  Catalog  of  Algorithms  and  Subroutines 
[1.1.2];  and  Part  IV;  User  Reference  Guide  [1.1.3]. 

CCFM. VENTS  adopts  and  consolidates  many  previously  established  zone-type  fire 
modeling  ideas.  As  the  discussion  proceeds,  this  report  will  place  the 
resulting  features  of  the  model  into  the  perspective  of  the  relevant  reference 
sources . 

1.2.  Generic  Features  of  the  CCFM 

Although  the  CCFM  would  be  developed  in  stages,  each  stage  would  have  certain 
common  generic  features.  These  are  identified  below. 


1.2.1.  The  Feature  of  Consolidation  of  Past/Future  Advances  in  Compartment 
Fire  Modeling  i 

The  feature  of  consolidation  means  that  during  the  CCFM  development  process 
full  advantage  would  be  taken  of  the  significant  number  of  proven  advances  of 
fire  modeling  research.  While  improvements  and/or  generalizations  of 


^Numbers  in  brackets  always  refer  to  the  list  of  references  at  the  end  of 
this  report. 
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available  products  may  be  required,  there  should  be  no  unnecessary  duplication 
of  fundamentally  sound  and  proven  technology. 


1.2.2.  The  Feature  of  Being  Well -Documented 

Quality  documentation,  which  would  be  critical  to  user-friendliness  of  CCFM 
products,  must  provide;  details  of  the  mathematical  algorithms/submodels  and 
associated  modular  subroutines  for  the  various  physical  phenomena  simulated; 
details  of  the  numerics  module (s)  used  to  solve  the  governing  set  of  the  model 
equations;  information  on  the  structure  of  the  overall  code  and  the  details  of 
any  utility  software  subroutines  used  to  manage  code  simulations;  and  a well- 
written  user  guide  which  complements  easy-to-use  input  and  output  software. 


1.2.3.  The  Feature  of  Modularity:  Code  Can  Receive  Advanced/Improved 
Submodels,  Solver (s ) , and  Input/Output  Software 

The  feature  of  modularity  would  be  critical  to  insuring  the  required  orderly 
and  efficient  future  growth  of  and  improvement  to  existing  stages  of  CCFM 
products.  The  use  of  well-defined  modular  elements  in  the  CCFM  would  also  be 
critical  for  a successful  team  development  effort.  Indeed,  by  constructing 
CCFM  products  with  modular  components  having  well-defined  inputs  and  outputs, 
a development- team  effort  that  even  involves  significant  inter- institution 
participation  would  be  possible. 


1.2.4.  The  Feature  of  User-Friendliness 

User-friendliness  refers  to  the  ability  of  people  with  a reasonably  wide  range 
of  technical  background,  other  than  the  members  of  the  CCFM  development  team, 
to  have  ready  access  to  CCFM  products.  Access  here  refers  to  the  ability  of  a 
serious  user  to:  run  any  particular  CCFM  code  application;  modify  the  code  as 
required;  and  extract  and  use  modular  portions  of  the  code  for  the  user's  oxm 
particular  needs.  User-friendliness  requires  easy-to-use  software  which  can 
be  run  on  generally  accessible  computer  hardware  and  which  is  supported  by 
quality  documentation. 


1.2.5.  The  Feature  of  Flexibility  - From  Basic  to  "Reference"  Modeling 

Successively  developed  CCFM  "applications"  would  lead  to  a capability  for 
simulating  a particular  fire  scenario  with  more  or  fewer  fire  - generated 
phenomena  and  with  different  levels  and  preferences  of  modeling 
sophistication.  Such  a range  of  capability  would  result  in  different  kinds  of 
user- flexibility . For  example,  a simulation  which  accounts  for  the  complete 
array  of  available  phenomena  and  which  is  treated  at  the  highest  possible 
level  of  modeling  sophistication  would  be  relatively  costly,  perhaps 
unnecessarily  so,  in  terms  of  required  user  and  computer  resources.  If  a user 
does  not  require  the  maximum  detail  and  sophistication,  or  is  not  prepared  to 
pay  the  costs,  then  simpler  or  more  basic  CCFM  simulations  would  be  available. 
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1.2.6.  The  Feature  of  Inter -Institution/International  Participation  and 
Co-operation  in  the  Development  of  the  CCFM 

There  would  be  clear  benefits  to  significant  inter- institution  and 
international  co-operation  in  the  development  of  the  CCFM.  These  would 
include;  the  pooling  of  resources,  the  avoidance  of  unnecessary  duplication  of 
effort,  and  the  generation  of  consensus  CCFM  products  which  would  maximize 
appropriate  proliferation  and  use.  It  is  envisaged  that  such  co-operation 
would  be  implemented  by  a variety  of  mechanisms  which  are  now  being 
identified. 


1.3.  Strategy  for  Achieving  a Final  Stage  of  CCFM  Development 

The  strategy  for  developing  the  CCFM  to  its  final  stage  includes  the 
following  components: 

• Develop  the  CCFM  in  stages. 

• Each  CCFM  stage  would  culminate  in  a progressively  more 
sophisticated,  well - documented , multi-room  fire  model  computer 
code . 

• Staged  products  would  involve  progressive  consolidation  of 
existing  and  developing  fire  modeling  technology.  As  required, 
appropriate  further  development  and  improvement  of  this  technology 
would  be  carried  out. 

• Development  of  a prototype/first  stage  CCFM  product:  CCFM. VENTS 
(i.e.,  the  Consolidated  Compartment  Fire  Model  application 
designated  as  VENTS). 

• Development  of  a second  stage  CCFM  product  (possible  name 
CCFM . BUILDING)  with  significant  advances  over  CCFM. VENTS  in 
modeling  of  combustion,  flow  dynamics,  and  heat  transfer;  not  yet 
of  "Reference  quality." 

• Anticipate  final,  "Reference"  level  of  sophistication  in  CCFM  by 
third- to- fourth  stage. 

1.4.  CCFM. VENTS:  The  Prototype/First  Stage  Application  of  the  CCFM  Concept 

In  terms  of  the  above  strategy,  this  documentation  and  the  associated  computer 
software  for  CCFM. VENTS  completes  the  prototype/first  stage  of  CCFM 
development . 

The  basic  features  of  CCFM. VENTS  include: 

• multiple  rooms/vents, 

• concentrations  of  0^^  and  other  products  of  combustion. 
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high  or  low  cross -vent  pressure  difference, 


forced  ventilation  via  simple  fan-duct  systems  - user  can  specify 
flow  rate  or  fan/duct  characteristics, 

multiple,  specified,  time  - dependent  fires, 

simplest  plume/heat  transfer, 

stack  effect, 

wind  loading, 

robust  solver, 

menu- driven  user-friendly  input, 

operation  on  computers  that  support  FORTRAN  77,  e.g.,  IBM- 
compatible  and  Apple  Mac  II  personal  computers, 

complete  modularity  of  all  physical -phenomena  algorithms  and 
associated  subroutines, 

complete  descriptions  of  all  physical  submodel  algorithms 
including  a catalog  of  all  algorithms  and  associated  subroutines, 

complete  description  of  the  structure  of  the  code  and  of  all 
associated  computer  software, 

user  guide. 


CCFM. VENTS  uses  a generic  CCFM  equation  set  that  would  be  applicable  for  all 
stages  of  the  CCFM  development  process.  This  equation  set  will  be  presented 
in  Section  2. 

A concise  but  complete  description  of  the  physical  basis  of  all  of  the 
physical  modeling  features  included  in  CCFM. VENTS  is  presented  in  Section  3. 
The  material  in  this  section  constitutes  the  basis  of  the  complete  CCFM.VENTS- 
specific  equation  set. 

Section  4 describes  the  method  of  specifying  a problem  for  a particular 
CCFM. VENTS  simulation.  As  will  be  seen,  the  varied  means  of  specifying  a 
problem  within  the  context  of  the  CCFM. VENTS  inputs  lead  to  a capability  for 
the  simulation  of  the  phenomena  of  wind  loading  and  stack  effect. 

This  Part  I presentation  concludes  with  Section  6,  which  provides  a summary 
and  overview  of  Parts  II -IV  of  this  report. 
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2.  The  Generic  CCFM  Equation  Set 


2.1.  Model  Assumptions  and  Some  Definitions 

The  CCFM  uses  a two-layer  zone-type  description  of  the  fire-generated 
environment  for  all  rooms  of  a facility.  This  means  that  the  environment  in 
each  room  is  assumed  to  be  comprised  of  one  or  two  relatively- quiescent 
stratified  zones  which  together  fill  essentially  all  of  the  room's  fixed 
volume.  These  zones,  depicted  in  Figure  2.1,  include  an  elevated  temperature 
buoyant  upper  layer  and  a relatively  cool  lower  layer. 

It  is  of  interest  to  note  that  there  are  a wide  range  of  previously  developed 
compartment  fire  models  which  use  a two -layer  zone -type  description  of  the 
fire  environment.  These  include:  CIT  [2.1.1],  Harvard  V [2.1.2],  RFIRES 
[2.1.3],  ASET  [2.1.4],  Harvard  VI  [2.1.5],  BRI  [2.1.6],  FAST  [2.1.7], 
NBS/Harvard  VI  [2.1.8],  DSLAYV  [2.1.9],  FIRST  [2.1.10],  CIFI  [2.1.11],  and 
MRFC  [2.1.12].  While  extensive,  the  latter  list  is  not  intended  to  be 
exhaustive . 

Both  layers  of  a room  are  assumed  to  be  spatially  uniform  in  temperature, 
density,  and  composition.  During  the  course  of  a simulated  fire  scenario, 
conservation  of  mass  and  chemical  species  are  imposed  on  all  layers  of  all 
rooms  of  a modeled  facility.  By  imposing  conservation  laws  and  specifying  an 
equation  of  state  for  the  layers,  assumed  here  to  be  consistent  with  a perfect 
gas  model  for  air,  one  is  led  to  a generic  set  of  model  equations  for  the 
instantaneous  layer  properties  in  an  arbitrary  room.  The  CCFM  generic  set  is 
derived  in  Appendix  A and  presented  in  Section  2.3.  Except  for  variations  in 
nomenclature,  choices  of  solution  variables,  and  questions  related  to 
treatment  of  pressure -work  terms  in  the  energy  equation,  there  must  be  a 
fundamental  equivalence  between  the  generic  layer-equation  sets  of  different 
zone -models . 

The  properties  of  the  layers  change  as  a function  of  time  by  virtue  of  mass, 
energy,  and  species  transfers  from  relatively  small-volume  plumes,  jets,  near- 
boundary flows,  combustion  zones,  and  other  isolated  or  distributed  sources. 
The  volumes  of  all  of  these  small-volume  zones  are  assumed  to  be  negligible. 
The  number  and  types  of  phenomena  occurring  in  these  zones  are  critical  to  the 
overall  model.  Indeed,  the  number  of  physical  and  chemical  phenomena  taken 
into  account,  and  the  detail  and  sophistication  with  which  such  phenomena  are 
modeled  actually  distinguish  one  overall  zone  fire  model  from  another  or,  in 
the  present  case,  one  application  of  the  CCFM  generic  equation  set  from 
another.  The  model  equations  for  all  of  the  source  and  transfer  zones  are 
designated  as  the  submodels  of  the  overall  model.  They  must  of  course  all  be 
internally  consistent  with  conservation  laws. 

Without  specifying  the  number  or  details  of  the  submodels  of  a particular 
version  or  application  of  the  CCFM,  the  following  definitions  and  designations 
are  useful  generally  for  deriving  and  implementing  the  generic  CCFM  equation 
set  for  a particular  CCFM  application: 

i i]  u i t L i 1 rates  of  mass  and  product  of 

combustion  k flowing  to  the  upper  [lower]  layer  of  room  i,  respectively. 
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Also,  Qu  ^ [Ql  i]  is  the  net  rate  of  enthalpy  plus  heat  transfer  plus  energy 
release  flowing  to  the  upper  [lower]  layer  of  room  i. 

Note  that  all  of  the  above  net  flow  rates  represent  the  sum  of  transfers  to 
the  layers  from  all  above-mentioned  plumes,  jets,  near-boundary  flows, 
combustion  zones,  and  other  isolated  or  distributed  sources  considered  and 
taken  into  account  in  any  particular  CCFM  application. 


2.2.  Pressure  in  the  Generic  Layer  Equations 

The  generic  CCFM  equations  will  take  into  account  pressure -work  terms  in  the 
energy  equation.  In  this  regard,  this  section  will  clarify  the  use  of 
approximate  representations  of  the  pressure  variable. 

Define  the  dimensionless  quantity 


5 = pg(y  - y)/p(y) 


(2.2.1) 


where  y is  an  arbitrary  elevation  within  a room  of  interest  above  a datum 
elevation,  y is  an  arbitrary  fixed  value  of  y within  the  room,  p is  the 
density  of  the  layer  at  y,  g is  the  acceleration  of  gravity,  and  p(y)  is 
absolute  pressure.  In  fire  scenarios  it  is  always  true  that 

|5|  « 1 (2.2.2) 

The  two-layer  modeling  assumption  and  Eq.  (2.2.2)  along  with  the  momentum 
equation,  i.e.,  hydrostatics,  and  the  perfect  gas  law  lead  to 

P(y)  = P(y)[l  + 0(5)]  ~ p(y)  (2.2.3) 


Note  that  although  pressure  differences  within  a room  are  negligible  the 
pressure  differences  between  one  room  and  another  may  not  be  negligible. 

The  Eq.  (2.2.3)  approximation,  adopted  implicitly  in  all  two-layer  zone 
models,  is  also  adopted  here.  In  every  room  i,  y is  taken  to  be  the  elevation 
of  the  floor,  ypLooR  i P(y)  is  designated  as  Pfloor  i-  Thus, 


Pfloor,!  ^ P(y  = yFLOOR.i)  ^ (2.2.4) 

Eqs . (2.2.3)  and  (2.2.4)  are  used  to  estimate  pressure-work  terms  in  the 
energy  equation  to  a consistent  level  of  accuracy  when  such  terms  are  non- 
negligible . 
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As  will  be  discussed  in  Section  3.3,  calculation  of  p(y)  to  determine  flow 
through  vents  due  to  cross -vent  pressure  differences  will  usually  require 
accurate  estimates  of  the  0(5)  term  in  Eq . (2.2.3).  Therefore,  when 
calculating  vent  flows  the  approximation  p(y)  = p(y)  of  Eq • (2.2.3),  used  to 
describe  thermodynamic  state  variables,  will  not  be  applicable. 


2.3.  Generic  Equation  Set  for  An  Overall  Facility 
2.3.1.  The  Generic  Equation  Set  for  An  Arbitrary  Room 

Consider  room  i of  a simulated  facility.  Define  as  the  fixed  volume  of  the 
room  and  ^ as  the  instantaneous  volume  of  its  upper  layer,  i.e.,  ^ 

is  the  volume  of  the  lower  layer.  Also,  consistent  with  the  definitions  in 
Section  2.1,  My  ^ , M^  ^ and  u i ’ L i defined  as  the  total 

instantaneous  mass  and  amount-of -product  k in  the  upper  and  lower  layers, 
respectively . 

Typical  fire  scenarios  involve  PpLooR  i which  vary  relatively  little 

from  a datum  absolute  pressure,  Pqatum  • example,  if  Poatum  ^ problem 

is  taken  as  the  local  atmospheric  pressure  at  ground  elevation,  of  the  order 
of  10^  pascals,  then  typical  variations  of  the  Pfloor  from  Pd^tum 

only  be  of  the  order  of  10^  pascals  or  less.  For  this  reason,  use  of  the 
transformed  pressure  variable, 


'^Pfloor, i Pfloor,!  ‘ Pdatum  (2.3.1) 

will  lead  typically  to  a higher  level  of  accuracy  in  the  computation  of 
Pfloor  i-  perturbation-pressure,  ^Pfloor  i’  adopted  as  a solution 

variable  in  place  of  the  pressure  variable,  Pfloor  i-  [Note  that  the  S of  Eq. 
(2.3.1)  is  part  of  the  single  symbol  5p,  and  should  not  be  confused  with  the  5 
of  Eqs.  (2.2.1)-(2.2.3) . ] 

In  the  Appendix,  the  generic  CCFM  equation  set  for  an  arbitrary  room  i has 
been  derived  from  application  of  the  laws  of  conservation  of  mass,  energy  and 
species.  As  mentioned,  pressure -work  terms  are  taken  account  of  in  the  energy 
equation.  The  equations  are  presented  in  terms  of  the  solution  variables 

, i ’ \,i’  , U , i ’ ^k  , L , i > , i ’ ^PfLOOR,!’ 

For  rooms  i of  a facility  with  horizontal  cross-sectional  areas,  A^ , which  are 
uniform  with  elevation,  (e.g.,  for  right  cylindrical  rooms),  the  solution 
variable  Vy  ^ can  be  replaced  by  the  alternate  solution  variable  yLAYER  i > 
elevation  of  the  interface  separating  the  upper  and  lower  layers. 

The  equations  from  the  Appendix  are  reproduced  below. 
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Mass  in  Upper  and  Lower  Lavers: 


dMu  ^/dt  = Mu  , (2.3.2) 

dMu_i/dt  = ML_i  (2.3.3) 


where,  consistent  with  the  definitions  in  Section  2.1,  the  terms  M^  ^ and  M^ 


individual  mass  flow  rate  contributions  to  the  layers  from  all  plumes,  jets, 
sources,  etc.  The  terms  u i ’ l i ’ Qu  i > i > to  be  introduced 

below,  will  represent  analogous  product  and  enthalpy  flow  rate  contributions, 
respectively . 

Volume  of  Upper  Laver  or  Elevation  of  the  Interface  (For  Rooms  With  Vertical 

Walls) : 


dVu_,/dt  = 

{(7-l)/[7(5PKL00R,i  +Pdatum)])[(1  - Vo,i/V,)Qu,i  - (V^ , , /V, )Q, , , 

(2.3.4) 


where 


7 = Cp/C^  = 1 + R/C,  (2.3.5) 

and  where  and  are  the  specific  heats  at  constant  pressure  and  volume, 
respectively,  R is  the  gas  constant.  All  of  these  properties  are 
approximateed  by  those  of  air. 

For  rooms  with  vertical  walls,  Eq . (2.3.4)  can  be  replaced  by 

*^yLAYER  , i ^ /(VPfloOR  , i [ (Yce  I L I NG  , i " Yl  A Y E R , i ^ , i 

^YlaYER,!  ' YfLOOR , i , i ^ 
(2.3.4') 


where  Yceiling  i Yfloor  i fixed  elevations  of  the  ceiling  and 

floor  of  room  i,  respectively. 
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Perturbation  Pressure  at  the  Floor : 


L OOR  , i Kt*  1)/^!  ] (Qu  , i Ql  , i ^ 

Amount  of  Product  k in  Upper  and  Lower  Lavers: 

dPk,u,i/dt  = Pk.u.i 
, L , i ^ , L , i 

Equation  of  State : 

^PfLOOR.I  PdATUM  ""  Pu,i'^u,i^  ^ 


(2.3.6) 


(2.3.7) 

(2.3.8) 


(2.3.9) 


Layer  Densities.  Temperatures,  and  Concentrations  of  Products: 


Py  ^ and  pj^  ^ , the  densities  of  the  upper  and  lower  layer,  respectively  are 
computed  according  to  their  definitions  from  the  solution  variables  My  ^ , 


Ml , i . and  Vy  _ . 


Pu.i  = Mu.i/Vy,, 

Pl.x  = Ml.l/CV,  - Vy_  . ) 


(2.3.10) 

(2.3.11) 


where 


(^CEILING , i ' yFLOOR.i^ 

, i ^ (yCE I L I NG , i * yLAYER.i^ 


(2.3.12) 

(2.3.13) 


With  the  above  values  of  the  densities,  the  absolute  temperatures  of  the  upper 
and  lower  layer,  Ty  ^ and  Ty  respectively  are  then  computed  from  the 
solution  variable  Spp^QQj^  ^ using  Eq . (2.3.9). 
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FLOOR , i 


(2.3.14) 


T 


U , i 


(5p 


Pdatum  (Pu  , 


'^L  , i (^PpLOOR.i  Pdatum  , i 


(2.3.15) 


Finally,  ^ ^ , Cj^  l i - concentrations  of  product  of  combustion  k in  the 
upper  and  lower  layer,  respectively  are  computed  from  the  solution  variables 
according  to  their  definitions. 


•k  , U , i ^k  , U , i /^U  , i 

(2.3.16) 

•k  , L , i ^ ^k  , L , i /^L  , i 

(2.3.17) 

The  final  set  of  generic  equations  for  the  solution  variables  is  Eqs . (2.3.2)- 
(2.3.8)  . Evaluation  of  the  terms  My  ^ , M^  ^ , Pj^  y ^ , Pj^  ^ ^ , Qy  ^ , and  Qy  y 
which  appear  on  the  right  sides  of  these  equations  will  be  discussed  in  the 
next  section.  The  derived  variables,  Py  ^ , Py  y , Ty  y , Ty  y , Cj^  y y , and 
L i ’ computed  from  the  solution  variables  by  using  Eqs.  (2.3.10)- 

(2’.  3’.17)  . 

2.3.2.  Forming  the  Equation  Set  for  the  Facility 

To  use  the  CCFM  to  simulate  a fire  scenario  it  is  generally  necessary  to  solve 
simultaneously  the  complete  equation  set  for  all  solution  variables,  described 
in  Section  2.3.1,  for  each  room,  i,  of  the  facility.  In  such  an  equation  set, 
the  equations  for  the  solution  variables  My  i > i ’ ^u  i ytAYER  i ’ 

rooms  i with  vertical  walls),  5ppyQQj^  y will  always  be  coupled  (i.e.,  the 
variables  appear  implicitly  or  explicitly  in  equations  other  than  their  own) , 
while  the  equations  for  the  u i ’ ^k  L i ’ k = 1 to  Npj^Qy  may  be 

uncoupled  (i.e.,  the  variables  do  not  appear  implicitly  or  explicitly  in  any 
equation  other  than  their  own).  Thus,  if  the  number  of  rooms  in  a simulated 
facility  is  , the  number  of  coupled  solution  equations  in  the  overall 

differential  equation  set  will  be  between  ^‘Nroom  ^prod^’^oom- 

the  present  context,  O2  would  be  considered  as  one  of  the  products  of 

combustion.)  As  will  be  seen,  in  the  present  CCFM  application  code, 

CCFM. VENTS,  the  differential  equations  for  the  P^  y y > and  P^  l i always 

uncoupled  and  the  number  of  coupled  differential  equations  is  4-Nj^qq(^. 

For  CCFM. VENTS,  the  generic  equations  are  implemented  with  the  assumption  that 
all  rooms  of  a modeled  facility  are  right  cylinders,  i.e.,  Eqs.  (2.3.12)  and 
(2.3.13),  with  elevation- independent  Ay,  are  applicable  and  Eq . (2.3.4')  is 
used  instead  of  Eq . (2.3.4). 

The  generic  equation  set  for  the  fire  environment  describes  explicitly  the 
state  of  the  two  gas/smoke  layers  in  each  room  of  the  facility.  Heat  and  mass 
transfer  exchanges  between  the  surfaces  of  the  solid  partition  elements  of  the 
rooms  and  the  gas  layers  are  included  implicitly  in  the  terms  My  y , My  y , 

^k  u i ’ ^k  L i > Qu  i ’ ^L  i '"^hich  appear  on  the  right-hand  sides  of  the 

various  components  of  this  equation  set.  In  any  particular  CCFM  application. 
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a treatment  of  these  terms  would  include  any  required  analysis  of  the 
diffusion  of  heat  and  mass  through  solids.  In  CCFM. VENTS,  heaf  and  mass 
diffusion  phenomena  are  not  treated  explicitly  and  such  analysis  is  not 
required.  To  explicitly  include  these  diffusion  phenomena  an  enriched 
structure  of  the  generic  equation  set  could  be  appropriate. 


2.4.  The  Equation  Set  for  a Specific  CCFM  Application 

One  establishes  a specific  application  of  the  generic  CCFM  equations  by 
identifying  explicitly  the  algorithms  used  to  calculate  the  terms  My  ^ , M^  ^ , 
u i ’ L i ’ Qu  i ’ i ^hich  appear  on  their  right-hand  sides.  As 

mentioned  in  Section  2.1,  the  actual  phenomena  included  in  a simulation  and 
the  detail  and  sophistication  with  which  they  are  modeled  would  distinguish 
one  application  of  the  CCFM  generic  equation  set  from  another.  The  remainder 
of  this  work  provides  the  details  of  the  first-stage  CCFM  application, 

CCFM. VENTS . 

The  next  section  identifies  and  presents  the  algorithms  for  all  the  phenomena 
taken  into  account  in  CCFM. VENTS. 
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3.  The  CCFM. VENTS -Specific  Equation  Set:  Submodel  Algorithms 

3.1.  The  Layers'  Sources  of  Mass,  Enthalpy,  and  Products  of  Combustion 

In  CCFM. VENTS  it  is  assumed  that  there  are  three  basic  phenomena  which  can 
lead  to  contributions  to  the  rates  of  flow  of  mass,  enthalpy  and  products  of 
combustion  to  the  upper  and  lower  layers  of  the  rooms  of  a facility,  namely, 
(1)  a fire  and  its  plume;  (2)  natural  flows  through  vertical  vents  driven  by 
room-to-room  cross-vent  hydrostatic  pressure  differences;  and  (3)  forced- 
ventilation  flows  between  pairs  of  rooms  connected  by  fan/duct  systems . 
Algorithms  used  to  model  these  phenomena  in  CCFM. VENTS  are  described  in  detail 
in  Sections  3.2,  3.3,  and  3.4,  respectively. 


3.2.  Plume  and  Heat  Transfer  in  Rooms  of  Fire  Origin 
3.2.1.  Some  Basic  Assumptions 

Consider  one  of  the  rooms  of  the  facility  which  contains  a fire.  This  is 
illustrated  in  Figure  3.2.1.  Figure  3.2.2  depicts  the  fire  and  its  plume  for 
all  possible  configurations  during  the  two- layer  simulation  of  the  fire 
environment.  Section  3.2  describes  an  algorithm  for  estimating  the  rates  of 
exchange  and  addition  of  mass,  enthalpy,  oxygen,  and  other  products  of 
combustion  between  and  to  the  two  layers.  The  exchanges  considered  here  are 
those  which  are  the  result  of  direct  action  of  the  fire/plume  dynamics.  The 
algorithm  is  based  on  the  model  of  reference  [3.2.1]  for  the  upward  mass  flow 
rate  in  a buoyant  point-source-driven  plume.  It  is  the  same  as  the  algorithm 
used  in  the  one-room  fire  model,  ASET  [2.1.4].  Thus,  it  is  assumed  that  the 
energy  and  all  products  of  combustion  released  by  the  fire  are  deposited  in 
the  room  by  means  of  a small  volume  source.  This  is  modeled  as  a point 
source . 

The  total  energy- release  rate  of  the  fire,  Q,  is  specified.  It  is  assumed 
that  the  combustion  zone  and  the  plume  radiate  from  the  source  at  an  effective 
rate  Aj^Q,  where  is  specified;  thus 


Qsource  ^ 


(3.2.1) 


is  the  part  of  the  fire's  energy- release  - rate  which  drives  the  convective 
buoyant  fire  plume  upward. 

Flaming  fires  exhibit  Aj, ' s in  the  range  0 < Ap,  < 0.6,  e.g.,  smaller  values  for 
smaller  methane  fires  and  higher  values  for  large  polystyrene  fires.  However, 
for  hazardous  fires  over  a wide  range  of  energy  release  rate  and  involving  a 
wide  range  of  common  groupings  of  combustibles,  it  is  reasonable  to 
approximate  flame  radiation  by  choosing  Aj^  = 0.35  [2.1.4].  This  is  taken  as 
the  default  value  in  CCFM. VENTS. 

The  fire  is  assumed  to  be  a source  of  oxygen  and  of  (NpjjQp  - 1)  other  user- 
specified  products  of  combustion  which  are  taken  account  of  in  the  simulation. 
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Following  Che  rule  of  reference  [3.2.2],  which  is  useful  for  most  fuels 
commonly  encountered  in  fires,  it  is  assumed  that  the  rate  of  O2  generation  in 
the  oxidation  process,  Pq2  source’  negative  (i.e.,  oxygen  is  consumed)  and 
directly  proportional  to  Q. 


rate  of  O2  generation  = P02, source 


SOURCE 


-F-Q 


(3.2.2) 


where 


F = 0.076(10-^ ) (kg  of  02)/J 


(3.2.3) 


It  is  assumed  that  the  other  products  of  combustion  are  generated  at  specified 
rates  source*  ^ = 2 to  Npj^Qp  [(unit  of  product  k)/s].  Typically,  the  rate 
at  which  fuel  mass  is  introduced  to  the  combustion  process  is  negligible 
compared  to  other  characteristic  mass  flow  rates  of  interest  in  the 
simulation,  and  this  is  assumed  to  be  the  case  here. 

If  the  lower  layer  is  of  non-zero  thickness  and  if  the  fire  is  in  the  lower 
layer,  then  it  is  assumed  that  a complete  description  of  the  lower- layer  room 
environment  at  the  instant  of  interest  is  known.  This  includes  the  density, 
absolute  temperature,  concentration  of  O2  [(kg  of  02)/(kg  layer)],  and 
concentrations  of  the  other  (Np^g^  - 1)  products  of  combustion  [ (unit  of 
product  k)/(kg  layer)].  These  are  designated  by  p-^  , T^  , Cq2  l ^ L’ 

Cj^  L > ^ 2 to  NppQp  , respectively.  The  O2  and  other  product  concentrations 

would  be  required  only  when  these  concentrations  were  actually  being 
predicted. 

The  ideal  gas  relationship  for  air  is  assumed  to  be  valid  for  the  entire  room 
environment . 

The  elevation  of  the  interface  separating  the  upper  and  lower  layers  is 
ytAYER  • fire  is  located  at  elevation  ypiRg,  which  is  assumed  to  be 

between  or  at  the  elevation  of  the  room's  floor  and  ceiling,  ypLOOR  ycEiL’ 

respectively.  Consistent  with  the  idealized  point-source  plume  model  of 
[3.2.1],  it  is  appropriate  to  choose  ypjRj;  as  the  lowest  elevation  above  the 
base  of  the  combustion  zone  where  relatively  free,  unrestricted,  lateral 
entrainment  of  the  local  gas  into  the  combustion  zone  or  plume  is  possible. 

Note  that  all  elevations  are  relative  to  the  same  datum  elevation. 

It  is  assumed  that  the  total  fraction  of  Q which  is  lost  to  the  bounding 
surfaces  of  the  room  by  all  modes  of  heat  transfer  is  given  by  A.j,Q,  where 
is  specified.  The  remaining  portion  of  Q is  Qxotal ’ 


0 T A L ' '^x)Q 


(3.2.4) 
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This  is  assumed  to  be  deposited  entirely  in  the  upper  layer  (except  for  the 
special  Configurations  3 and  6,  Figure  3.2.2,  to  be  discussed  below).  When 
the  fire  is  in  the  upper  layer  this  is  all  added  to  the  upper  layer  directly 
and  there  is  no  other  enthalpy  flow  due  to  the  fire.  When  the  fire  is  in  the 
lower  layer,  the  net  rate  of  enthalpy  flow  into  the  upper  layer  due  to  the 
fire  is  Qxotal  rate  of  enthalpy  flow  associated  with  the  gases 

entrained  into  the  plume  from  the  lower  layer  and  deposited  by  the  plume  into 
the  upper  layer 

Some  limited  guidance  on  the  choice  of  in  the  range  0.6  < A^  < 0.9  is 
provided  in  reference [ 2 . 1 . 4 ] . The  default  value  of  A^  is  taken  to  be  0.8  in 
CCFM. VENTS . 

The  algorithm  assumes  no  gas  - to  - room- surface  losses  of  products  of  combustion. 

The  fire -driven  exchange  and  addition  of  mass,  enthalpy,  and  products  of 
combustion  between  and  to  the  two  layers  can  depend  on  the  strength  of  the 
fire  plume  and  on  the  relative  elevations  of  fire  and  interface  and  interface 
and  ceiling.  At  any  instant,  a fire  scenario  can  fall  into  any  one  of  the  six 
configurations  shown  in  Figure  3.2.2.  The  flows  for  each  of  these 
configurations  are  discussed  below. 

3.2.2.  Fire  in  the  Lower  Layer 

3. 2. 2.1.  Interface  Below  the  Ceiling 

Configuration  1: 

This  configuration  satisfies 


(3.2.5) 


ytAYER  ■ ypiRE  Z > 0 and  ycEir  ^ yiAYER 


and  is  depicted  in  Figure  3.2.2a. 

Following  reference  [3.2.1],  the  mass  flow  rate  in  the  plume  as  it  penetrates 
the  layer  interface  and  enters  the  upper  layer  is  estimated  to  be 


(3.2.6) 


(3.2.7) 


where  Cp  is  the  specific  heat  of  air  at  constant  pressure,  assumed  to  be 
constant,  and  g is  the  acceleration  of  gravity. 


15 


MpLUME  is  modeled  as  being  identical  to  the  flow  extracted  by  plume 
entrainment  from  the  lower  layer.  Accordingly,  the  fire  plume  leads  to  the 
following  contributions  to  the  total  rates  of  flow  to  the  lower-layer  of  mass, 
; enthalpy,  : molecular  oxygen,  l ^ i*!  L’  other  products  of 

combustion,  Pj^  ^ ^ = 2 to  Npj^Qp  : 


M = - M 
^ ^ L u M E 


Ql  ■ ^PLUME  ■ *2p  • 


P = - c • M 

■^02  , L ^02  , L ^^PLUME 


P = - C -M 

, L k , L “plume 


, k = 2 to  N 


PROD 


(3.2.8) 

(3.2.9) 

(3.2.10) 

(3.2.11) 


At  the  interface,  all  of  the  mass,  enthalpy,  and  products  of  combustion 
flowing  in  the  plume,  both  from  generation  in  the  fire  and  from  plume 
entrainment,  are  deposited  into  the  upper  layer.  Accordingly,  the  fire  and 
its  plume  leads  to  the  following  contributions  to  the  total  rates  of  flow  to 
the  upper-layer  of  mass.  My;  enthalpy,  ; molecular  oxygen,  P^^  u ^ u> 
other  products  of  combustion,  Pj^  ^ ; 


Mu 

= M 

“plume 

(3.2.12) 

Qu 

“ MpmME‘^p''^L  Qtotal 

(3.2.13) 

^0  2 

.U  ‘^0  2,L‘MpLUME  ■ 

(3.2.14) 

Pk. 

U ~ ^k.L'MpLUME  ^k,  SOURCE’  ^ ~ ^ ^pROD 

(3.2.15) 

3. 2. 2. 2.  Interface  at  Ceiling,  Zero-Thickness  Upper  Layer 

This  configuration  is  associated  with  an  instant  of  time  during  the  fire 
scenario  when  the  upper  layer  has  not  yet  formed  or  when  a previously- 
developed  upper  layer  has  been  reduced  to  one  of  zero  thickness.  At  such  a 
time  it  is  important  to  determine,  within  the  context  of  the  overall  two -layer 
zone -type  approximation,  whether  the  characteristic  temperature  of  the  plume 
at  the  ceiling  above  the  lower- layer  temperature  is  large  enough  to  initiate 
the  growth  of  a coherent  stratified  upper  layer.  If  so,  then  it  would  be 
reasonable  to  carry  out  all  flow- exchange  calculations  as  in  Configuration  1. 
This  would  lead  to  the  initiation  of  upper- layer  growth.  If  the 
characteristic  temperature  rise  of  the  plume  is  smaller  than  some  specified 
value,  ATgpg , then  it  is  reasonable  to  assume  that  all  of  the  plume  flow  is 
mixed  back  into  the  full - room-height  lower  layer.  Under  such  a circumstance, 
no  mass,  enthalpy,  or  products  in  the  plume  flow  would  be  available  to 
initiate  growth  of  an  upper  layer. 
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An  appropriate  specification  of  AT^-pg  would  depend  on  the  fire  scenario  being 
simulated.  For  example,  pre-ignition  stratification  of  a room,  by  whatever 
means  (e.g.,  due  to  solar  heating  effects),  can  establish  near-ceiling 
temperatures  which  are  several  degrees  K above  the  average  room  temperature. 
With  such  an  initial  environment  it  is  inappropriate  to  suggest  that  the 
environment  as  modified  by  the  fire  can  be  described  usefully  by  a two- layer 
model,  unless  the  difference  between  the  predicted  upper  and  lower  layer 
temperatures  is  greater  than  the  characteristic  temperature  difference  of  the 
pre- ignition  stratified  environment.  For  example,  if  the  pre- ignition 
temperature  difference  is  estimated  to  be  5-10  K,  then  it  may  be  reasonable  to 
choose  ATppg  as,  e.g.,  10-20  K. 

In  general,  it  is  inappropriate  to  take  seriously  a two- layer  description  of 
room  fire  environment  unless  the  difference  in  temperature  between  the  layers 
is  at  least  of  the  order,  of  1 K. 

The  ATgpg-type  of  temperature  - rise  criterion  for  initiation  of  an  upper  layer 
is  adopted  here,  where  the  value  of  AT^pg  in  CCFM. VENTS  is  taken  to  be  1 K. 

The  characteristic  plume  temperature,  Tp^uME > taken  to  be  the  average 
temperature  of  the  plume  gases  at  the  ceiling  elevation,  and  the 
characteristic  temperature  rise,  ATp^ui^^g  , used  in  the  criterion  is  found  from 


where 


AT 

“■^PLUME 

^SOURCE 


I'^PLUME  ■ QsOURCe/^*^p‘^PLUMe) 

and  MpLUME  calculated  in  Eqs . (3.2.1), 


(3.2.16) 

(3.2.6)  and  (3.2.7). 


Configuration  2:  Relatively  Strong  Plume  Initiates  Upper  Layer  Growth 
This  configuration  satisfies 

Flayer  ‘ ^fire  = Z > 0 and  = Ylayer  (3.2.17) 

^'^PLUME  — ^"^EPS  (3.2.18) 


and  is  depicted  in  Figure  3.2.2b. 

The  condition  of  Eq . (3.2.18)  indicates  that  the  plume  at  the  ceiling  is 
relatively  strong.  It  is  therefore  assumed  that  growth  of  an  upper  layer  will 
be  initiated.  Under  these  circumstances,  all  flow- exchange  calculations  are 
carried  out  as  in  Eqs.  ( 3 . 2 . 8 ) - ( 3 . 2 . 15 ) of  Configuration  1. 
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Configuration  3:  Relatively  Weak  Plume  Does  Not  Initiate  Upper  Layer  Growth 


This  configuration  satisfies 

ytAYER  ■ ypiRE  = Z > 0 and  y^giL  ^ ytAYER 

(3.2.19) 

^"^PLUME  '^'^EPS 

(3.2.20) 

and  is  depicted  in  Figure  3.2.2c. 

The  condition  of  Eq . (3.2.20)  indicates  that  the  plume  at  the  ceiling  is 
relatively  weak.  It  is  therefore  assumed  that  growth  of  an  upper  layer  is  not 
initiated  as  a result  of  fire-generated  flows.  Under  these  circumstances,  the 


fire  simply  deposits  its  products  into  the  lower  layer, 
following  contributions  to  the  rates  of  change  of  mass , 
in  the  upper  and  lower  layers 

This  leads  to  the 
enthalpy,  and  products 

^ ^0  2,U  ^ , U ^ Qu  = 0>  k = 2 to  Npj^Qp 

(3.2.21) 

Ql  ^ Qtotal 

(3.2.22) 

P02,L  = - F-Q 

(3.2.23) 

P =P  k=2toN 

^k  , L -^k,  SOURCE’  lu  i-Np^Qp 

(3.2.24) 

3.2.3.  Fire  in  Upper  Layer 

Configuration  4: 

This  configuration  satisfies 

ytAYER  ■ ypiRE  = Z < 0 

(3.2.25) 

and  is  depicted  in  Figure  3. 2. 2d. 

As  indicated  in  the  figure,  although  a fire  plume  will  be  formed  in  the  upper 
layer  it  is  assumed  that  the  plume  does  not  lead  to  any  exchanges  of  mass  and 
products  between  the  layers.  The  fire  simply  deposits  its  products  into  the 
upper  layer.  This  leads  to  the  following  contributions  to  the  rates  of  change 
of  mass,  enthalpy,  and  products  in  the  upper  and  lower  layers 
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(3.2.26) 


(3.2.27) 


P, 


02  , U 


F-Q 


(3.2.28) 


P, 


k , U 


P, 


k , SOURCE  ’ 


k = 2 to  N- 


PROD 


(3.2.29) 


3.2.4.  Fire  at  the  Layer  Interface 

During  the  course  of  a simulation  it  is  possible  that  the  elevations  of  the 
fire  and  layer  interface  coincide  exactly.  When  this  occurs,  yL^YER  ' ^fire 
Z = 0 and,  according  to  Eqs . (3.2.6)  and  (3.2.7),  MpLUME  ^ subcases 

this  configuration  have  to  be  considered. 

3. 2. 4.1.  Fire  Below  the  Ceiling 

Configuration  5: 


This  configuration  satisfies 


(3.2.30) 


Flayer  ' ^fire  ^ ^ yFiRE  ycEiLiNc 


and  is  depicted  in  Figure  3 . 2 . 2e . 

This  is  the  case  where  there  is  a non- zero  upper  layer  and  the  elevation  of 
the  fire  and  of  the  interface  are  identical. 

From  Eqs.  (3.2.6)  and  (3.2.7)  fipmME  ^ Then  the  flow  contributions 
predicted  by  the  equations  of  Configurations  1 [Eqs.  (3.2.8)-(3.2.15)]  and  ■ 
[Eqs.  ( 3 . 2 . 26 ) - ( 3 . 2 . 29 ) ] are  identical  and  either  set  is  applicable. 

3. 2. 4. 2.  Fire  at  the  Ceiling 

Configuration  6: 

This  configuration  satisfies 


(3.2.31) 
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and  is  depicted  in  Figure  3.2.2f. 

This  is  the  case  where  there  is  a zero  upper  layer  and  the  elevations  of  the 
fire,  the  interface,  and  the  ceiling  are  all  identical.  Under  these 
circumstances  and  for  the  plume  model  adopted  here,  no  plume  exists  to  entrain 
mass  from  the  lower  layer  and  to  initiate  growth  of  an  upper  layer.  It  is 
therefore  assumed  that,  as  in  Configuration  3,  growth  of  an  upper  layer  will 
not  be  initiated  as  a result  of  fire  - generated  flows  and  that  the  fire  simply 
deposits  its  products  into  the  lower  layer.  This  leads  to  contributions  to 
the  rates  of  change  of  mass,  enthalpy,  and  products  in  the  upper  and  lower 
layers  which  are  identical  to  those  prescribed  in  Eqs . ( 3 . 2 . 21) - ( 3 . 2 . 24) . 


3.2.5.  The  Algorithm  PLUGO  and  Its  Associated  Subroutine 

The  algorithm  described  above  in  Sections  3.2.1  - 3.2.4  for  calculating  fire- 
plume  - generated  flows  and  heat  transfer  has  been  named  PLUGO.  A concise 
presentation  of  PLUGO  and  its  associated  FORTRAN  subroutine  appears  in  Part 
III  [1.1.2]  of  this  work. 

3.3.  Natural  Vent  Flows  and  Heat  Transfer 
3.3.1.  Three  Stages  of  the  Calculation 

Consider  a pair  of  adjacent  rooms  of  a facility  involved  a fire  scenario. 

Room  pairs  considered  here  are  those  which  have  a common  vertical  wall  segment 
with  an  arbitrary,  non- zero  number  of  penetrations,  referred  to  here  as  vents. 
Described  in  Section  3.3  is  an  algorithm  for  estimating  the  rates  of  exchange 
of  mass,  enthalpy,  oxygen,  and  other  products  of  combustion  through  the  vents 
and  between  the  layers  of  the  two  rooms.  The  exchanges  considered  here  are 
the  result  of  vent  flows  driven  only  by  room-to-room  cross-vent  hydrostatic 
pressure  differences.  Only  rectangular  vents  with  horizontal  and  vertical 
boundaries  are  considered. 

At  any  instant,  the  complete  vent- flow  calculation  is  carried  out  in  three 
stages.  In  the  first  stage,  all  parameters  are  established  for  a calculation 
of  the  room-to-room  hydrostatic  pressure  difference  at  an  arbitrary  elevation 
of  the  common-wall  segment.  This  pressure  difference  drives  the  room-to-room 
vent  flows  to  be  studied  here.  In  the  second  and  third  stages  of  the 
calculation,  the  individual  vents  in  the  wall  segment  are  considered  one  by 
one.  The  second  stage  determines  all  characteristics  of  any  particular  vent 
flow.  As  will  be  discussed  below,  such  a flow  can  involve  from  one  to  six, 
contiguous,  horizontal  slabs  of  uniform-property,  unidirectional  flows.  For 
each  such  slab,  the  third  stage  determines  the  portion  of  the  flow  going  to 
each  of.  the  two  receiving-room  layers  and  the  rate  of  extraction  from  its 
single  source-room  layer. 

The  natural  vent  flow  and  heat  transfer  model  described  in  Section  3.3  is 
based  on  a consolidation  of  and  extensions  to  ideas  presented  in  references 
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[2.1.1],  [2.1.2],  [2.1.4],  and  [ 3 . 3 . 1 - 3 . 3 . 5 ] . These  references  will  be 

identified  at  appropriate  places  in  the  text.  It  should  also  be  noted  that 
the  basic  ideas  used  here  for  calculating  vent  flow  rates  under  typical 
conditions  of  relatively  low  cross-vent  pressure  difference  have  also  been 
used  in  all  vent  flow  calculations  of  the  models  of  references  [2.1.1-2.1.12]. 


3.3.2.  Specification  of  Geometry  and  Fire  Environment  of  the  Room  Pair 

Designate  two  adjacent  rooms  of  interest  as  rooms  1 and  2,  and  refer  to  Figure 
3.3.1.  The  index  i will  be  used  below  to  indicate  physical  characteristics 
and  layer  properties  of  room  i,  i = 1 or  2 . 

The  geometric  characteristics  of  the  two  rooms  required  for  the  present 
calculation  are  the  elevations  of  their  floors  and  ceilings,  Yfloor  i 
ycEiL,i>  respectively. 

In  each  room  the  instantaneous  states  of  the  layers  are  specified  by  the 
elevation  of  the  interfaces,  Ylayer  i>  separating  the  upper  and  lower  layers 
and  by  the  properties  of  the  layers’  ^ , Pl  , i . , i - . i > ^02  . u , i ’ ^02  , l , i ’ 

u i ’ ^k  L i ’ ^ ^ ^ ^PROD  ’ defined  in  Section  2. 

At  instants  of  time  when  there  is  no  upper  layer  [lower  layer],  i.e.,  when 
YLAYER.i  = ycEiL.i  [yLAYER.i  =yFLOOR.i]-  Upper-layer  [lower-layer]  properties 
are  not  used  in  the  calculation.  Also,  as  in  Section  3.2,  O2  and  other 
product  concentrations  would  be  required  only  when  the  simulation  under 
consideration  is  actually  being  used  to  predict  these  concentrations. 

The  final  specifications  of  the  environment  are  the  ^Pp^ooR  i 
rooms.  As  defined  in  Eq . (2.3.1),  these  are  the  differences  between  the 
absolute  hydrostatic  pressures  at  the  floor  of  each  room  and  some  fixed, 
absolute,  hydrostatic  datum  pressure,  Pdaxum • discussed  in  Section  2.3.1, 

the  5ppLooR  i solution  pressure  variables  used  in  the  basic  equation 

set  for  room  i. 

One  of  the  two  rooms  connected  by  a vent  can  be  an  outside -environment  space 
or  "room."  Local  to  any  single  outside  - exposed  vent,  the  hydrostatic  pressure 
in  the  outside  space  is  computed  in  the  same  way  as  in  the  inside  room.  In 
the  outside  space,  however,  the  vent  is  always  submerged  by  a "lower  layer"  of 
unlimited  extent  and  of  uniform,  specified  properties  analogous  to  p^  ^ , T^  ^ , 
C02  L i > ^k  L i ’ . the  layer- interface  elevation,  Ylayer  i ’ 

the  outside  space  or  "room"  is  indefinitely  large.  In  order  to  fix  the 
pressure  calculation  in  the  outside  space,  a reference  outside-space 
elevation,  analogous  to  Yfloor  i the  inside  rooms,  and  a pressure  at  this 

reference  elevation,  analogous  to  5pFL00R,i  for  the  inside  rooms,  must  be 
specified. 

In  the  discussion  to  follow,  when  one  of  the  facilitY  "rooms"  is  an  outside 
space,  room  variables  identified  should  be  substituted  bY  analogous  variables 
of  the  outside  space. 
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As  in  Section  3.2,  in  Section  3.3  all  elevations  are  relative  to  the  same 
datum  elevation. 


3.3.3.  First  Stage  of  the  Calculation:  Parameters  to  Determine  the  Cross-Vent 
Pressure  Difference  at  an  Arbitrary  Elevation  of  the  Common-wall  Segment 

At  an  arbitrary  elevation,  y,  in  one  of  the  rooms,  (y)  is  calculated  from 
Pi  (y)  = ^Pfloor,!  Pdatum  (3.3.1) 


where  5p^ (y) , the  difference  between  the  absolute  hydrostatic  pressure  at  y 
and  at  ypLooR  i ’ obtained  from  hydrostatics  to  be 


Pt.iS^y  ■ yFLOOR,i)>  yFLOOR.i  - y - yLAYER,i> 


5p.  (y)  = ^ 


Pl , i S^ytAYER , i 


(3.3.2) 

yFLOOR.i^  ■ Pu,iS(y  " yLAYER.i^’ 
yLAYER.i  y - yCEIL.i 


and  where  g is  the  acceleration  of  gravity. 

In  general,  the  above  distributions,  5p^ (y) , would  involve  two  linear 
functions  of  elevation  for  each  of  the  two  rooms.  This  is  depicted  at  the  top 
of  Figure  3.3.2.  One  function  describes  8p^(y)  in  the  lower  layer,  from  the 
floor  of  the  room  to  the  interface,  and  the  other  in  the  upper  layer,  from  the 
interface  to  the  ceiling.  At  the  interface,  5p^ (y)  is  continuous,  but  its 
slope  is  discontinuous. 

At  a given  y of  the  common-wall  segment,  the  difference  between  the  pressures 
in  the  two  rooms,  Apj^  2>  is  defined  by 

APi,2  ^ Pi(y)  - P2(y)  (3.3.3) 

Using  Eqs  . (2.3.1),  ( 3 . 3 . 1)  - ( 3 . 3 . 3 ) , it  is  found  that  Ap^^  ^ ^^.n  be  calculated 
from 

^Pi  , 2 ^ (^Pfloor,!  ' ^PFLooR.a^  [^Pi(y)  ' '5p2(y)]>  (3.3.4) 

ywiN  ^^^(yFLOOR  , 1 > yFL00R,2^  - y - ™i’^(yCEIL  , 1 ’ yCEIL,2)  ^ yMAX 
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Note  that  Eq . (3.3.4)  is  only  valid  at  elevations  at  or  between  the  lower-  and 
upper-elevation  limits  of  the  common-wall  segment,  and  y[.5y^x  > 

respectively . 

The  Eq . ( 3 . 3 . 4)  - formulation  for  calculating  Apj^  ^ reveals  the  benefit  of  using 
the  relative  pressure,  Spp^QQj^  instead  of  the  absolute  pressure,  Pj^  (y  = 
ypLOOR  i ) ’ solution  pressure  variable  in  the  basic  equation  set 

discussed  in  Section  2.  If  p^^  (y  = ypLooR  used,  then  Ap^^  ^ would  have 

to  be  calculated  directly  from  Eq . (3.3.3).  This  would  Involve  the  difference 
of  the  two  absolute  pressures,  p^  (y  = ypLOOR  Pa  (y  ^ ypLOOR  • Since 

it  would  be  typical  for  these  two  values  to  be  identical  to  each  other  to 
within  four  or  five  significant  figures,  there  would  be  a corresponding  loss 
in  the  number  of  significant  figures  in  the  computed  value  of  Ap^  i.e.,  as 
mentioned  earlier,  absolute  pressures  are  typically  of  the  order  of  100,000 
pascals  while  pressure  differences  throughout  entire  facilities  are  typically 
of  the  order  of  a few  or,  at  most,  several  tens  of  pascals.  This  would,  in 
turn,  lead  to  an  unnecessary  loss  in  accuracy  in  the  vent- flow  rates,  to  be 
determined  below  from  Ap^  2 • By  using  ^Pfloor  i solution  pressure 

variable  and  by  introducing  an  appropriate  value  for  Pq^jum  iu  Eq.  (3.3.1), 
the  significant  figures  which  would  be  lost  when  calculating  Ap^  ^ from  Eq . 
(3.3.3)  is  regained  by  calculating  it  from  Eq . (3.3.4). 

The  bracketed  term  of  Eq . (3.3.4)  involves  the  difference  of  the  two  5pj^  (y) 
functions  where  each  one  involves  up  to  two  linear  segments  continuous  in  the 
interval  y[.5iN  < y < y^^x  • Therefore,  Apj^  ^ is  also  continuous  in  this 
interval  and  is  made  up  of  from  one  to  three  linear  segments.  The  vertical 
bounds  of  these  segments  involve  a minimum  of  two  and  a maximum  of  four 
discrete  elevations:  at  the  bottom,  y^Ax  ’ top.  and,  if  they  lie  in 

between,  either  or  both  of  the  yLAYER  i • It  is  a simple  matter  to  calculate 
from  Eq . (3.3.4)  the  value  of  Ap^^  ^ which  corresponds  to  each  of  these  N^^J.LEV 
wall  elevations,  where  2 < N^j^ev  — Assuming  that  this  is  done,  the  unique 

set  of  points,  [y^  . Ap^^  2 ^yN  ^ 1 > N = 1 to  , which  define  the  Ap^  ^ 

segments  is  then  ordered  according  to  monotonically  increasing  yj^  . 

In  general,  each  of  the  one  to  three  linear  segments  of  Ap^^  ^ can  have  a slope 
of  arbitrary  magnitude  and  sign.  For  this  reason,  it  is  possible  for  each  to 
have  a zero  Ap^^  ^ value  inside  its  vertical  y-range  of  applicability.  For  any 
particular  segment,  it  is  a simple  matter  to  determine  if  such  a zero  exists. 
This  could  be  done,  for  example,  by  establishing  whether  or  not  the  now-known 
Apj^  2 values  at  the  segment  endpoints  are  both  non-zero  and  of  different  sign. 
If  it  is  determined  that  a segment  has  an  elevation  where  Ap^^  ^ =0,  then  this 
elevation  is  calculated  and,  together  with  its  corresponding  Ap^^  2 ^ value, 

it  is  sorted  into  the  previous  [y^  . Ap^^  2 ^yN  ^ ] set.  After  all  zeros  have  been 
determined,  the  resulting  [y^,  , Ap^^  2 ^yN  ^ ^ would  involve  points, 

where  N^j-^ev  satisfies  2 < N^j-^ev  - 

The  final  N^g^Ev  points  identifies  end-points  of  one  to  six,  uniform- 

sign,  contiguous  segments  of  Ap^^  ^ . These  segments  are  linear  in  y and  they 
determine  completely  the  distribution  of  Ap^^  ^ which  is  continuous  across  the 
entire  common-wall  segment. 
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Figure  3.3.2  illustrates  possible  ^ distributions  with  from  one  to  six 

segments  and  with  corresponding  endpoints  , Ap^^  2 ^yN  ^ ^ > 1 ^ N < N^elev  ^ 
to  7 . 


3.3.4.  Second  Stage  of  the  Calculation:  Characteristics  of  a Vent  Flow 
3. 3. 4.1.  The  Uniform- Property  Flow  Slabs 

Consider  an  arbitrary  vent  in  the  above  common-wall  segment.  Its  width  and 
upper  and  lower  bounds  are  specified  as  b,  yvxop  yvBOX  > respectively , 

where  yvjQp  - yvBoi  ^ cross-vent  Ap^^  2 distribution  is  that  portion  of 

the  Ap^  2 distribution  for  the  common-wall  segment,  given  by  Eq . (3.3.4), 
which  lies  in  the  interval  yvBox  - y — yvxop • 

To  obtain  the  end-points  of  the  Ap^^  ^ distribution  for  the  vent:  use  Eq . 
(3.3.4)  to  calculate  the  two  Apj^  2’^  corresponding  to  yvBox  yvxop  > sort 

the  resulting  two  [y,  Ap^  2]  points  into  the  [y^  , Ap^^  2 ^yw  ^ ^ 
entire  common-wall  segment,  and  eliminate  any  points  of  the  resulting  set 
having  < yysox  yw  - yvxop  • Non-unique  points  of  the  set  are  also 

eliminated,  e.g.,  if  yvBox  ^ ypLOOR  i-  resulting  set  of  vent  elevations 

would  have  points,  where  N^^lev  satisfies  2 < N^j-^ev  - 

The  final  N^g^EV  points  identify  end-points  of  one  to  six,  uniform- sign , 
contiguous  segments  of  Ap^  ^ which  determine  the  distribution  of  Ap^^  ^ across 
the  vent . 

At  a given  y,  the  direction  and  rate  at  which  flow  is  driven  through  the  vent 
is  determined  locally  by  the  sign  and  magnitude  of  Ap^^  ^ . If  Ap^^  ^ > Q>  for  a 
given  segment  of  the  Ap^  ^ distribution,  then,  throughout  this  segment,  the 
flow  will  be  from  room  1 to  room  2 and  the  source  of  the  flow  will  be  one  of 
the  one  or  two  uniform-property  layers  of  room  1.  If  Apj^  2 < ^ ^he 

segment,  then  the  flow  is  reversed  and  its  source  will  be  one  of  the  one  or 
two  layers  of  room  2.  The  uniform-property  flow  corresponding  to  a single 
Apj^  2 segment  is  referred  to  here  as  a flow  slab. 

For  a generic  vent.  Figure  3.3.3  illustrates  possible  Ap^^  ^ distributions  with 
from  one  to  six  segments  and  N^j-^ev  corresponding  endpoints.  In  the  figure, 
the  direction  of  the  flow- slab  velocity,  corresponding  to  each  Ap^  ^ segment, 
and  actual  velocity  distributions  are  indicated  by  arrows. 


3. 3. 4. 2.  Mass  Flow  Rate  in  the  Flow  Slab 

Number  the  flow  slabs  through  a vent  from  the  lowest  (N  = 1)  to  the  highest  (N 
^ ^SLAB  ^VELEV  ' 1^6).  Then,  according  to  the  above  - implemented  sorting 
procedure  the  bottom  and  top  elevation  of  an  arbitrary  slab,  N,  will  be  y^,  and 
respectively.  Attention  here  is  confined  to  the  problem  of  estimating 
the  rate  of  cross -vent  flow  in  slab  N. 

Between  y^^  and  y^  + i.  the  direction  of  the  flow  in  the  slab  and  the  uniform 
density,  temperature,  and  concentrations  of  products  of  combustion  of  the  flow 
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upstream  of  the  vent  are  fixed.  Also,  throughout  the  depth  of  the  slab,  p^ , 

P2  , and  Apj^  ^ each  linear  in  y.  Ap^^  ^ t)e  zero  at  the  bottom  or  top  of 

a slab,  and  it  can  be  identically  zero,  i.e.,  a slab  with  zero  cross-vent 
flow,  but  Ap]^  2 cannot  change  signs  in  the  slab. 

Identify  the  flow  direction  of  the  slab  by  =1,  0,  or  -1,  according  to 

whether  the  flow  is  from  room  1 to  room  2 (Apj^  2 > ^ in  the  slab)  , zero 
(Apj^  2 = bhe  slab)  , or  room  2 to  room  1 (Apj^  2 ^ in  the  slab)  , 

respectively.  If  the  flow  in  the  slab  is  zero  then  no  further  consideration 

is  required.  Therefore,  assume  that  the  arbitrary  slab  under  consideration 
has  a non-zero  flow. 

With  complete  generality,  assume  that  the  flow  in  the  slab  is  from  room  1 to 
room  2.  (If  the  flow  is  actually  from  room  2 to  room  1,  then  re -number  for 
now  the  pair  of  rooms  under  consideration.)  Also,  designate  the  properties  of 
the  slab's  room-1  source  layer  with  the  subscript  1,  i.e.,  the  density, 
temperature,  oxygen  and  product  k concentrations  are  p^,  , Cq2  1,  and  c^^ 

respectively . 

Define  Ap  and  Ap  according  to 


[\p  = p^  - p^  (3.3.5) 

Ap(y)  = Ap^_2(y)  = Pi (y)  - P2(y)  - ^ (3.3.6) 

where  p^  and  P2  are  the  uniform  densities  in  rooms  1 and  2,  respectively,  at 
the  height  of  the  slab.  (Recall  that  while  Eq . (3.3.6)  is  similar  to  Eq . 
(3.3.3),  the  rooms  have  been  re-numbered  so  that  the  inequality  of  Eq . (3.3.6) 
is  valid.)  Define  Apjj  , pj^  and  Ap^j^^^,  Pj^  + i as  the  previously  obtained  values 
of  Ap , Pj^  at  y = yj,  and  y^  + i.  respectively.  Define  also 


e(y)  = Ap(y)/p^(y)  = 1 - x(y)  > 0 

x(y)  ^ P2(y)/Pi(y)  = i - «(y) 


(3.3.7) 


In  typical  compartment  fire  scenarios,  cross -vent  pressure  differences,  which 
are  relatively  small  (compared  to  p^^ ) and  lead  to  relatively  low-speed  flow 
velocities  (compard  to  the  speed  of  sound) , satisfy 


0 < Ap/p^  = e (y)  « 1 for  "typical"  compartment  fires  (3.3.8) 


Fire  scenarios  involving  flow  slabs  with  relatively  large  cross -vent  pressures 
and  associated  high-speed  flow  velocities  will  have  e's  which  do  not  satisfy 
Eq.  (3.3.8).  Such  moderate-size  e's  will  lead  to  single-slab  vent  flows  of 
the  type  depicted  in  Figure  3.3.3(a).  The  Ap  distributions  for  such  moderate - 
e flows  will  satisfy  | Apj,  - Apj^^  ^ | /min(Ap^  , Apj,^^)  « 1. 
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Define  y as  the  elevation  of  the  middle  of  the  slab,  and  e as  the  value  of  e 
at  5'  , i . e . , 


y ^ (Yn  + yN+i)/2 

i = €(y  = y)  = 1 - X = (AP(^  + Apj,+  i)/(p^  + Pn+i) 


(3.3.9) 


Assuming  isentropic  flow  and  adopting  the  traditional  approximation  that  the 
thermodynamic  state  of  the  gas  flowing  through  the  vent  varies,  at  a given  y, 
only  in  the  streamwise  direction,  the  total  rate  of  mass  flow  in  the  slab, 
Mslab’  estimated  from  (see  reference  [3.3.1]) 

^SLAB  = bCw[  (8p^)i^V3]  (y^  + i - Yn)- 

[(I^PnI  + |APnAPn  + i I + |Ap„^  J )/(  |Apj^  I + |Ap^^ji/2)j 

(3.3.10) 


where 


w = w(x)  = w(l  - e)  ; Z = C(x)  = C(1  - i) 


(3.3.11) 


and  where  w and  C are  given  by 


w(x)  = f(x)/[2(l  - x)]^/2 


(3.3.12) 


f(x) 


{[27/(7  - 1)]x2/t'[1  - X > [2/(7  + 

(7[2/(7  + i)](7+i)/(7-i) }l/2  X < [2/(7  + 1)]7/(7-i) 

(3.3.13) 


C(x)  = 0.68x  + 0.85(1  - x)  = 0.85  - 0.17x  ; 0 < x < 1 (3.3.14) 


In  Eq . (3.3.13),  7,  the  ratio  of  specific  heats  of  the  flowing  gas,  taken  to 
have  the  properties  of  air,  is  1.40. 

In  the  above,  C is  a flow  coefficient  which  depends  generally  on  Reynolds 
number  and  x.  For  perfectly  isentropic  flow,  C = 1.  For  real  nozzles,  C can 
be  very  close  to  1 , of  the  order  of  0.99,  for  large  "pipe"  Reynolds  number 
flows  of  10®  or  more  [3.3.2,  p.99].  For  sharp-edged  orifices,  C = C(x)  and  is 
approximately  linear  with  x,  ranging  from  0.68  at  x = 1 to  0.85  at  x = 0 . 
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From  Eq . (3.3.7)  it  is  clear  that  a scenario  where  the  pressure  ratio  x is 

close  to  0 corresponds  to  a situation  where  the  absolute  pressure  in  room  1 is 
several  times  greater  than  that  in  room  2.  This  is  well  below  the  value  x = 
0.528,  required  to  initiate  choked  flow  through  the  vent.  For  orifice -like 
openings  typical  of  vents  in  compartment  fires,  x is  close  to  1 and  C = 0.68 
[3.3.3].  It  is  reasonable  to  approximate  C(x)  by  Eq . (3.3.14)  unless  a more 
precise  description  of  C is  known  for  some  particular  vent  configuration  of 
interest . 

The  following  expansion  is  useful  when  computing  w for  small  e 


lim  w(x)  = lim  w(l  - e)  = 1 - [3/(47)]e  + O(e^)  (3.3.15) 

x-^1  e-^0 


When  € 1 - X 0,  i.e.,  C ncompress  i ble  ^ Eq . (3.3.10)  is 

identical  to  the  small-cross-vent-pressure-difference  flow  equation  presented 
in  references  [3.3.3]  and  [2.1.2].  In  this  form  it  is  also  equivalent  to  the 
flow  calculations  described  in  [3.3.4]  and  in  most  of  the  models  of  [2.1.1]- 
[2.1.12].  Thus,  a simple  interpretation  of  Eq . (3.3.10)  is  that  when 
compressibility  effects  are  arbitrarily  strong  or  weak,  the  compartment  fire- 
generated flow  rate  of  a slab  of  flow,  originating  from  a uniform  stagnation 
source  and  passing  through  a constant-width  vertical  vent,  can  be  computed  by 
ignoring  elevation- related  differences  of  compressibility  in  the  slab  and  by 
using  the  characteristic  mass  flux  for  compressible  nozzle/orifice  flow  which 
corresponds  to  the  mid-slab  cross-vent  pressure  difference. 


3. 3. 4. 3.  Other  Flow  Rates  of  the  Flow  Slab 

With  Mslab  Eq . (3.3.10),  and  with  the  known  values  of  T^  , Cq2  i-  and 

Cj^  the  rates  of  flow  of  enthalpy,  oxygen,  and  product  of  combustion  k in 
the  slab,  which  are  designated  by  Qslab  > ^02  slab  > slab’  respectively, 

are  obtained  from 


Qslab  ^slab  ‘ ' '^1 

P = fi  • c 

^02  , SLAB  “slab  '-02,1 

P — M • c 

, SLAB  ^^SLAB  ^k  , 1 


(3.3.16) 

(3.3.17) 

(3.3.18) 


These  flows  come  from  the  source  layer  in  room  1.  Subsequent  to  the  slab  flow 
entry  to  room  2,  these  flows  are  divided  and  deposited  into  the  two  layers  in 
room  2 according  to  the  rule  to  be  specified  below  in  Section  3.3.5. 
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3. 3. 4. 4.  Properties  of  the  Penetrating  Flow 

In  order  to  use  Eq . (3.3.10)  in  any  particular  room  fire  model  it  is  necessary 

to  know  the  properties  of  the  penetrating  flow  as  it  enters  room  2.  With  such 
information  in  hand  it  would  then  be  possible  to  invoke  the  model's  adopted 
rule  for  dividing  and  depositing  this  incoming  flow  into  the  room's  upper  and 
lower  layers.  The  desired  result  is  obtained  by  observing  that  as  the  flow 
passes  through  the  penetrations  it  maintains  an  approximately  constant  total 
enthalpy,  i.e.,  the  flow  is  essentially  adiabatic.  Thus,  assuming 
deceleration  of  the  flow  to  a negligible  horizontal  velocity  prior  to  mixing 
in  the  receiving  room,  the  temperature  at  this  point,  Tpj.jj  , will  be  identical 


T 


PEN 


Ti 


(3.3.19) 


Using  this  result,  the  density  after  penetration  but  before  mixing, 
be  obtained  from  the  perfect  gas  law. 

^PEN  = Pi-(P2/Pi)  = (3.3.20) 

For  typical  fire  scenarios,  where  the  cross-vent  pressure  ratio,  x,  is  close 
to  1,  the  above  leads  to  the  expected  result  that  the  properties  of  the 
penetrating  flow  are  identical  to  its  properties  in  room  1. 

Conservation  of  species  leads  to 


^02  , PEN  ‘^02,1 
, PEN  ""  ^k  , 1 


(3.3.21) 


where,  again,  the  subscript  PEN  indicates  a property  of  the  slab  flow  in  the 
receiving  room  after  penetration  and  deceleration,  but  before  significant 
mixing . 


3. 3. 4. 5.  A Characteristic  Elevation  of  the  Flow  Slab 

As  the  flow  slab  enters  the  receiving  room,  possible  rules  for  its 
distribution  will  depend  generally  on  a characteristic  elevation  of  the  slab. 
Several  choices  for  this  elevation  are  available  directly,  viz.,  the  top, 
bottom,  or  midpoint.  Another  possible  choice  is  ysL^p , the  centroid  of  the 
distribution  of  slab  N's  entering  horizontal -velocity , or  -momentum  profile, 
i.e.. 
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(3.3.22) 


ySLAB 


^N+i  yN+1  yN+1  ^N+i 

= J pvydy  / J pvdy  = J vydy  / J vdy 

Yn  Yn  Yn  Yn 


where  p is  the  near-uniform  density  and  v = v(y)  is  the  velocity  distribution 
of  the  flow  slab  as  it  enters  room  2.  To  the  extent  that  the  horizontal 
penetration  of  a slab  into  the  receiving  room  is  dependent  on  its  horizontal 
momentum,  Yslab  ^^Y  t>e  preferable  to  other  characteristic  elevations  as  a 
measure  of  the  elevation  of  its  flow  as  it  is  projected  into  the  receiving 
room  by  its  horizontal  momentum. 

Throughout  the  height  of  the  slab 

v(y)  = A|Ap(y) 1^/2  (3.3.23) 


for  some  constant  A,  where 


Ap(y)  = A -I-  By  (3.3.24) 

In  Eq.  (3.3.24),  A and  B can  be  determined  from  the  known  values  of  yj^  , Apj^ 
and  y^i  + i.  ^Pn+1-  Then  y can  then  be  expressed  as  a function  of  Ap : 


Y = [Yn+i(^P  - ^Pn)  • Yn(^P  - ^Pn+1 ) ]/(^Pn+i  ' ^Pn ) (3.3.25) 

dy  = [(yN  + i ■ Yn)/(^Pn+i  ■ Apj,)]d(Ap)  (3.3.26) 

Eqs . (3.3.23),  (3.3.25),  and  (3.3.26)  are  now  used  to  change  the  variable  of 
integration  on  the  right  side  of  Eq.  (3.3.22)  from  y to  Ap . Carrying  out  the 
integration  leads  to 

ysLAB  - [(3/5)(Ap„,P'2  - = - y„)/(Ap„,P'"  - Ap„  = '2) 

+ (^PKtiyB  - *PByNtl)l/<'^PN-H  - “^Pn) 

(3.3.27) 


When  Apj^^j^  - Apj,  -*■  0,  i.e.,  the  velocity  profile  across  the  height  of  the  slab 
is  relatively  uniform,  an  analysis  of  Eq.  (3.3.27)  yields  the  expected  result 
that  the  centroid  of  the  momentum  profile  is  at  the  mid- elevation  of  the  slab, 
i.e.  , 
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YsLAB  (yN+l 

(^Pn  + 1 -^Pn  ) 0 


(3.3.28) 


IvThen  computing  Yslab  > • (3.3.28)  should  be  used  instead  of  Eq . (3.3.27)  when 

the  removable  singularity  of  Eq.  (3.3.27),  at  the  limit  (Apj^^  ^ -Apj^ ) 0 , is 

being  approached. 


3. 3. 4. 6.  Completing  the  Second  Stage  of  the  Calculation 

From  the  results  of  Sections  3 . 3 . 4 . 2 - 3 . 3 . 4 . 5 the  following  properties  of 
arbitrary  flow  slab  N,  can  be  calculated: 


^DIR’  ^SLAB’  QsLAB’  ^0  2,  SLAB’ 
, PEN 


k , SLAB 


ys 


LAB 


•PEN 


"PEN  ’ '-02  , PEN 


and 


The  second  stage  of  the  calculation  is  completed  by  calculating  this  array  of 
properties  for  each  of  the  slabs  of  the  vent  flow  under  consideration. 


3.3.5.  Third  Stage  of  the  Calculation:  Distribution  of  the  Slab  Flows  to  the 
Receiving  Layers  and  Rate  of  Extraction  from  the  Source  Layers 

3 . 3 . 5 . 1 . General  Aspects  of  the  Rule  for  Distributing  the  Slab  Flows 

The  third  and  final  stage  of  the  calculation  determines  the  contributions  of 
all  the  slab  flows  of  a vent  to  the  net  rates  of  flow  of  mass,  enthalpy, 
oxygen,  and  other  products  of  combustion  in  the  upper  and  lower  layer  of  each 
of  the  two  rooms.  These  are  the  net  rates  of  flow,  ^ , Qy  ^ , Pq2  y ^ , 

^k  u i ’ ^L  i ’ Ql  i > ^0  2 L i ’ ^k  L i ’ discussed  in  Section  2 and  required 

at  any  instant  in  the  calculation  of  the  right-hand-sides  of  Eqs . 2.3.2  - 
2.3.8  for  each  room  i of  the  simulation.  The  contributions  are  determined  one 
flow  slab  at  a time. 

A flow  slab  originates  from  a single  layer  of  the  source  room.  With  the  now- 
known  properties  of  the  flow  slab,  its  contributions  to  the  net  rates  of  flow 
in  the  source  room  can  immediately  be  taken  account  of.  For  example,  the 
instantaneous  net  rate  of  change  of  mass  of  the  source  layer  is  reduced  by 
Mslab  > while  that  of  the  other  layer  of  the  source  room  is  unaffected  by  the 
flow. 

As  a flow  enters  the  receiving  room  it  is  assumed  to  be  deposited  entirely 
into  one  of  the  two  layers  or  to  be  divided  between  them.  The  rule  for 
apportioning  the  flow  is  based  on  an  extension  of  a simple  scheme,  proposed  in 
reference  [3.3.5],  In  determining  the  net  enthalpy  flow  rate  deposition,  a 
rule  is  adopted  which  takes  account  of  heat  transfer  to  the  walls  of  the  room. 
This  rule  is  based  on  an  extension  of  ideas  proposed  in  reference  [2.1.1]  and 
used  in  reference  [2.1.4].  These  rules  will  be  described  below. 
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The  rule  of  apportioning  the  slab's  mass,  oxygen,  and  product  of  combustion 
flow  rates  is  based  on  the  temperature  of  the  penetrating  flow  relative  to  the 
temperature  of  the  receiving  room's  upper  and  lower  layers.  The  rule 
addresses  the  special  scenarios,  analogous  to  those  treated  in  Section 
3. 2. 2. 2,  when  it  is  not  reasonable  to  simulate  the  otherwise  expected 
initiation  of  upper-  or  lower- layer  growth  from  zero  thickness,  because  the 
magnitude  of  the  difference  between  the  slab  temperature  and  that  of  the 
existing  layer  is  smaller  than  some  specified  (relatively  small)  value,  AT^p^ . 

In  addition  to  the  latter  measures  of  relative  buoyancy,  the  rule  for 
apportioning  the  slab's  enthalpy  flow  rate  and  determining  the  heat  transfer 
to  the  surfaces  of  the  receiving  room  will  depend  explicitly  on  Tp , the 
temperature  of  the  layer  in  the  receiving  room  which  includes  the  elevation 

^SLAB  ’ Tp,  = Ty.iREC  ySLAB  ^ yLAYER.iREC’ 

Designate  the  fraction  of  the  flow  slab  deposited  into  the  upper  and  lower 
layers  as  and  , respectively.  Let  iREC  = 1 or  2 be  the  value  of  i for 
the  receiving  room. 

The  basic  rule  adopted  here  for  estimating  the  values  of  F^  and  F^  is 
described  below  for  three  possible  categories  of  interface  configuration. 


3. 3. 5. 2.  Interface  in  the  Receiving  Room  Between  the  Floor  and  the  Ceiling 

These  configurations  involve  upper  and  lower  layers  which  are  both  of  non- zero 
thickness.  They  satisfy 

yFLOOR.iREC  ytAYER.iREC  yCEIL.iREC  (3.3.29) 

Three  temperature  categories  are  distinguished  for  these  configurations. 

Temperature  Category  1:  Penetration  Temperature  Greater  than  that  of  the  Upper 

Layer 

This  category  satisfies: 


T 


PEN 


> T 


U , iREC 


(3.3.30) 


Here,  the  slab  flow  is  upward-buoyant  relative  to  the  upper  layer,  and  it  is 
assumed  that  the  entire  slab  flow,  with  its  oxygen  and  products  of  combustion 
is  deposited  there.  Thus, 


Fy  = 1 . and  F^^  = 0 . 


(3.3.31) 
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Temperature  Category  2:  Penetration  Temperature  Less  than  that  of  the  Lower 

Layer 


This  category  satisfies: 


"^PEN  ^ '^L.iREC 


(3.3.32) 


Here,  the  slab  flow  is  downward -buoyant  relative  to  the  lower  layer,  and  it  is 
assumed  that  its  entire  slab  flow,  with  its  oxygen  and  products  of  combustion 
is  deposited  there.  Thus, 


Fy  = 0 . and  = 1 . 


(3.3.33) 


Temperature  Category  3:  Penetration  Temperature  at  or  between  those  of  the 

Upper  and  Lower  Layers 


This  category  satisfies: 


T < T 

, iREC  - -^PEN 


< T 


U , iREC 


(3.3.34) 


Here,  the  slab  flow  is  neutral  or  upward -buoyant  relative  to  the  lower  layer, 
and  neutral  or  downward -buoyant  relative  to  the  upper  layer.  Here,  the 
adopted  apportionment  rule  is  that  the  slab  flow  is  divided  and  deposited  into 
the  two  layers  in  direct  proportion  to  the  temperature  differences.  In 
particular , 


Fu 

^"FpEN  ■ FL.iREC^/^Fu^iREc 

- Tl, 

iREC  ^ 

Fl 

('^U.iREC  ■ FpEN^/^'Fu.iREC 

- Tl, 

i REC  ) 

3. 3. 5. 3.  Zero -Thickness  Upper  Layer  in  the  Receiving  Room 

For  this  configuration  the  slab  elevation  is  always  within  the  elevation  range 
of  the  lower  layer  which  fills  the  room. 


ysL 


AB 


< 


AYER , iREC 


= yc 


EIL  , iREC 


(3.3.36) 
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Two  temperature  categories  are  distinguished. 


Temperature  Category  1:  Relatively  Strongly -Upward -Buoyant  Slab  Flow 


This  category  satisfies: 


0 < AT 


EPS 


< T 


PEN 


-L  , iREC 


(3.3.37) 


Here,  the  entire  slab  flow  is  deposited  in  a newly- initiated  upper  layer,  and 
Fy  and  F^  are  calculated  according  to  Eq . (3.3.31). 


Temperature  Category  2:  Relatively  Weakly-Upward- Buoyant , Neutrally-  or 

Downward- Buoyant  Slab  Flow 


This  category  satisfies: 


'^PEN  ■ "^L.iREC  '^'^E  P S 


(3.3.38) 


Here,  the  entire  slab  flow  is  deposited  into  the  existing  lower  layer  and  F^ 
and  Fj^  are  calculated  according  to  Eq . (3.3.33). 


3. 3. 5.4.  Zero -Thickness  Lower  Layer  in  the  Receiving  Room 

For  this  configuration  the  slab  elevation  is  always  within  the  elevation  range 
of  the  upper  layer  which  fills  the  room. 


Flayer, iREC  ypLOOR.iREC  ^ ysLAB  (3.3 

Two  temperature  categories  are  distinguished. 

Temperature  Catesorv  1:  Relatively  Strongly -Downward- Buoyant  Slab  Flow 


This  category  satisfies: 


0 < ATgpg  < - Tppjj 


(3.3.40) 
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Here,  the  entire  slab  flow  is  deposited  in  a newly- initiated  lower  layer,  and 
Fy  and  are  calculated  according  to  Eq . ( 3 . 3 . 33 ) . 


Temperature  Category  2:  Relatively  Weakly- Downward- Buoyant , Neutrally-  or 

Upward- Buoyant  Slab  Flow 


This  category  satisfies; 


'^U.iREC  ■ '^PEN  ^ ^"^EPS 


(3.3.41) 


Here,  the  entire  slab  flow  is  deposited  into  the  existing  upper  layer  and  F^ 
and  Fl  are  calculated  according  to  Eq . (3.3.31). 

3. 3. 5. 5.  Rule  for  Enthalpy  Deposition  and  Heat  Transfer  to  Room  Surfaces 

After  entering  the  receiving  room,  but  before  significant  mixing  there,  the 
slab  flow,  at  temperature  Tp^.^,  , is  surrounded  by  a layer  environment  at 
temperature  Tp . 

The  total  enthalpy  flow  rate  of  the  slab,  Qslab > first  assumed  to  be 
divided  and  distributed  to  the  upper  and  lower  layers  as  described  in  Sections 
3. 3. 5. 2 - 3. 3. 5. 4.  Then  it  is  assumed  that  a portion  of  Qsl^b  is  transferred 
to  the  surfaces  of  the  receiving  room.  The  portion  transferred  is  assumed  to 
be  a specified  fraction,  , of  the  "enthalpy  of  buoyancy,"  defined  as  the 
enthalpy  flow  rate  of  the  slab  calculated  relative  to  the  temperature 
difference  Tp^^,  - Tp  . 


"enthalpy  of  buoyancy"  = Mslab ' Cp • (TpEp  - Tp ) 

^Qslab'(^  ■ '^r/'^pen)  (3.3.42) 


This  "enthalpy  of  buoyancy"  of  the  slab  is  the  actual  source  of  relative  local 
upward-buoyancy  (if  the  enthalpy  of  buoyancy  is  > 0)  or  downward -buoyancy  (if 
the  enthalpy  of  buoyancy  is  < 0)  of  the  slab  in  the  receiving  room  at 
elevation  ysLAB • This  is  the  buoyancy  which  leads  to  the  distribution  of  the 
slab  flow,  by  convection,  according  to  the  rule  of  Sections  3. 3. 5. 2 - 3. 3. 5.4. 

The  "enthalpy  of  buoyancy"  would  tend  to  generate  an  upward-  or  downward- 
buoyant  plume  which  would  flow  toward  the  ceiling  or  floor,  respectively,  of 
the  receiving  room.  As  such  a plume  impinged  on  the  ceiling  or  floor  surface, 
it  would  lead  to  a ceiling  or  floor  jet.  This  would  lead  in  turn  to  a 
convective  heat  transfer  exchange  between  the  gas  and  the  surface.  The  heat 
transfer  phenomena  and  losses  for  an  upward-buoyant  enthalpy  of  buoyancy  is 
completely  analogous  to  the  heat  transfer  and  losses  for  a fire  plume.  For  a 
fire  plume,  the  enthalpy  converted  upward  from  the  fire's  combustion  zone,  (1 
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- \)Q  [see  Eq.  (3.2.1)],  is  directly  analogous  to  the  enthalpy  of  buoyancy  of 
the  slab.  Also,  for  an  idealized  fire  with  = 0,  i.e.,  with  negligible 
radiation  loss,  is  completely  analogous  to  [see  Eq . (3.2.4)]. 

It  is  consistent  with  the  guidance  on  the  choice  for  A^  in  reference  [2.1.4] 
to  choose  Xy  in  the  range  0.4  < Xy  < 0.8.  The  default  value  for  Xy  is  taken 
to  be  0.6  in  CCFM. VENTS. 

Following  the  above  ideas,  the  net  rate  of  heat  transfer  to  the  room  surfaces 
is 


net  rate  of  heat  transfer  to  surfaces  = -^v'Qslab'^^  ' ’^r/'^pen^ 

(3.3.43) 


and,  the  net  rate  of  enthalpy  flow  deposited  in  the  layers  of  the  receiving 
room  is 


net  rate  of  enthalpy  flow  to  layers  = Qslab  ' ^v'Qslab'^^  * "^r/'^pen^ 

^QsLAB't^  ■ ■ '^r/'^pen)] 

(3.3.44) 

The  rule  for  distributing  the  net  rate  of  enthalpy  flow  of  Eq . (3.3.44) 
between  the  two  layers  is  taken  to  correspond  to  the  earlier  rule  for 
distributing  mass,  oxygen,  and  products  of  combustion  according  to  the 
fractions  F^  and  F^  determined  in  Section  3. 3. 5. 2.  Thus, 


rate  of  enthalpy  flow  to  upper  layer  = Fu-QslabC^  ' ‘ '^r/'^pen^^ 

(3.3.45) 

rate  of  enthalpy  flow  to  lower  layer  = Fl ’Qslab  ‘ ’Fr/’^pen^] 

To  the  extent  that  an  estimate  of  the  actual  distribution  of  heat  transfer  to 
surfaces  is  required,  it  is  reasonable  to  assume  that  the  net  rate  of  heat 
transfer  to  room  surfaces  of  Eq . (3.3.43)  is  distributed  between  surfaces 

engulfed  by  the  upper  and  lower  layers  according  to  the  values  of  Fg  and  Fj^  . 
Thus , 
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rate  of  heat  transfer  to  upper  surfaces  = Fy  ' ’ Qs L ab  ^ ^ ' '^r/'^pen^ 

(3.3.46) 

rate  of  heat  transfer  to  lower  surfaces  = ’Qslab  " '^r/'^pen^ 

where,  again,  and  are  determined  as  in  Section  3. 3. 5. 2. 

It  is  noteworthy  that  the  above  estimates  lead  to  the  result  that  the  rates  of 
heat  transfer  to  surfaces  will  be  positive  or  negative  depending  on  whether 
Tpgj^  is  greater  or  less  than  , respectively.  As  illustrated  in  the  two 
examples  below,  this  is  consistent  with  the  assumption  that  in  some  semi- 
quantitative  sense  there  is  a correspondence  between  the  characteristic 
surface  temperature  of  the  room  and  . 

In  the  first  example,  assume  that  > '^r  > i-®.,  (1  - Tp/Tp^jj)  > 0,  and  that 

the  characteristic  temperature  of  the  upper  surfaces  of  the  receiving  room  is 
at  the  relatively- low- temperature , Tp  (e.g.,  as  in  the  situation  of  a room 
adjacent  to  a room  of  fire  origin  where  the  room-of -origin  upper  layer  is  just 
starting  to  be  transferred  through  the  upper  part  of  an  open  doorway) . Then 
the  slab  flow  will  lead  typically  to  a relatively  high- temperature  upward- 
buoyant  plume  in  the  receiving  room  which  will  impinge  on  the  relatively  low- 
temperature  ceiling.  As  predicted  in  Eqs . (3.3.45)  and  (3.3.46),  this  will 
lead  to  a positive  rate  of  heat  transfer  to  the  ceiling  (heating  of  the 
ceiling)  and  a corresponding  negative  component  of  rate  of  enthalpy  addition 
to  the  plume/upper  layer  (cooling  of  the  ceiling- impinging  plume). 

In  the  second  example,  assume  that  < Tp  , i.e.,  (1  - Tp/Tp^j^)  < 0,  and 

that  the  characteristic  temperature  of  the  lower  surfaces  of  the  receiving 
room  is  at  the  relatively  high- temperature , Tp  (e.g.,  as  in  a scenario  where 
relatively  cool  air  from  an  adjacent  space,  the  source  of  the  flow  slab,  is 
introduced  near  the  ceiling  of  the  receiving  room,  which  is  also  a room  of 
fire  origin) . Then  the  slab  flow  will  lead  typically  to  a relatively  low- 
temperature  downward -buoyant  plume  which  will  impinge  on  the  relatively  high- 
temperature  floor.  As  predicted  in  Eqs.  (3.3.45)  and  (3.3.46),  this  will  lead 
to  a negative  rate  of  heat  transfer  to  the  floor  (cooling  of  the  floor)  and  a 
corresponding  positive  component  of  rate  of  enthalpy  addition  to  the 
plume/upper  layer  (heating  of  the  floor- impinging  plume). 


3. 3. 5. 6.  Completing  the  Third  Stage  of  the  Calculation;  Natural -Vent -Flow 
Calculations  for  the  Total  Facility 

Contributions  of  any  slab  flow  to  the  upper  and  lower  layers  of  a pair  of 
vent  - connected  rooms  can  be  calculated  by  implementing  the  results  and 
discussion  of  Sections  3. 3. 5. 2 - 3. 3. 4. 5.  The  contribution  to  the  layers  from 
an  entire  single  vent  flow  is  then  obtained  by  carrying  out  sequentially  the 
slab  flow  calculations  for  all  of  the  one-to-six  slab  flows  of  the  vent.  This 
would  complete  the  third  stage  of  a vent- flow  calculation. 


36 


The  contribution  to  the  layers  from  all  vents  in  a single  common-wall  segment 
of  a room  pair  is  obtained  by  carrying  out  the  vent  calculations  sequentially 
for  all  vents  in  the  common-wall  segment.  This  would  involve  both  second-  and 
third- stage  types  of  calculations. 

The  contributions  of  natural  vertical  vent  flows  to  the  layers  of  all  rooms  of 
a facility  is  obtained  by  carrying  out  the  common-wall  segment  vent 
calculations  sequentially  for  all  vent- containing  common-wall  segments  of  the 
facility.  This  would  involve  all  three  stages  of  the  vent-flow  calculations 
of  Section  3.3. 

It  is  important  to  note  that  when  one  of  the  two  rooms  in  a calculation  is  the 
space  associated  with  the  outside  environment,  modifications  to  that  outside 
space  as  a result  of  leaving  or  entering  flow  slabs  are  ignored. 


3.3.6.  The  Algorithms  and  Associated  Subroutines,  COMWLl,  DELP,  VENTHP,  and 
FL0G02,  Used  to  Carry  Out  the  Three  Stages  of  the  Natural  Vent  Flow 
Calculation 

The  methodology  described  in  Sections  3.3.1  - 3.3.5  for  estimating  the  rates 
of  exchange  of  mass,  enthalpy,  oxygen,  and  other  products  of  combustion 
through  vertical  natural  vents  is  implemented  in  four  separate  algorithms 
named  DELP,  COMWLl,  VENTHP,  and  FL0G02 . 

DELP,  which  is  used  by  COMWLl  and  VENTHP,  is  a utility- type  of  algorithm  which 
implements  the  pressure  calculations  of  Eqs . (3.3.2)-(3.3.5).  COMWLl,  VENTHP, 
and  FL0G02  implement  the  calculations  of  Stage  1,  2,  and  3 and  their 
corresponding  Sections  3.3.3,  3.3.4,  and  3.3.5,  respectively.  The  use  of 
these  four  algorithms  is  depicted  in  Figure  3.3.4. 

A concise  presentation  of  DELP,  COMWLl,  VENTHP,  and  FL0G02  and  their 
associated  FORTRAN  subroutine  appears  in  Part  III  [1.1.2]  of  this  work. 


3.4.  Forced  Vent  Flows  and  Heat  Transfer  - Two  Options 

This  Section  describes  an  algorithm  for  estimating  the  rates  of  flow  of  mass, 
enthalpy,  oxygen,  and  other  products  of  combustion  through  a well - insulated , 
forced-ventilation , fan/duct  system  connecting  any  two  rooms,  or  a room  and 
the  outside  environment  in  a multi -room  facility.  The  rates  of  flow  to  and 
from  the  upper  and  lower  layers  of  two  connected  spaces  are  determined 
explicitly.  Figure  3.4.1  presents  a sketch  of  the  fire-generated  environment 
in  a pair  of  rooms  connected  by  such  a fan/duct  system. 

The  terminology  "well - insulated , " used  above  to  describe  the  fan/duct  system, 
is  meant  to  indicate  that  any  heat  transfer  to  the  walls  of  the  simulated 
fan/duct  systems  is  assumed  to  be  negligible. 

The  algorithm  is  implemented  with  one  of  two  basic  options.  In  Option  1,  the 
direction  of  flow  and  the  volume  flow  rate  through  the  system  is  specified. 

In  the  more  complicated  option.  Option  2,  the  instantaneous  flow  rate  and 
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direction  of  the  flow  are  determined  from  specified  fan  characteristics,  a 
specified  duct  resistance,  and  from  the  instantaneous  characteristics  of  the 
two- layer  room  environment  local  to  the  end-points  of  the  fan/duct  system 
being  simulated.  In  carrying  out  the  calculation  under  Option  2,  modeling  of 
the  fan/duct  ventilation  system  follows  an  extension  of  the  ideas  presented  in 
reference  [3.4.1], 


3.4.1.  Two  Stages  of  the  Calculation 

The  forced-vent  flow  calculation  is  carried  out  in  two  stages. 

In  the  first  stage,  the  rates  of  extraction  from  the  upper  and  lower  layers  of 
the  source  room  are  determined  along  with  all  the  properties  of  the  assumed- 
well  -mixed  the  duct  flow.  When  Option  2 is  implemented  and  system 
characteristics  are  used  to  determine  flow  direction  and  net  flow  rate,  an 
additional  algorithm  is  used  to  perform  the  actual  system-flow  simulation. 

This  algorithm  is  described  below  in  Section  3.4.5. 

The  first- stage  calculation  concludes  with  a description  of  the  flow  from  the 
source  room  into  the  duct,  in  the  duct  itself,  and  from  the  duct  into  the 
receiving  room.  The  flows  in  and  out  of  the  duct  are  described  in  the  manner 
of  the  VENTHP  natural -vent  flow  slabs  of  Section  3.3.4.  Here,  there  would  be 
only  one  or  two  entry  flow  slabs  and  one  or  two  exiting  flow  slabs  depending 
on  whether  or  not  a smoke  layer  interface  in  the  source  or  receiving  room  lies 
within  the  elevation  limits  of  the  fan/duct  system  vents  in  the  respective 
rooms . 

In  the  second-stage  calculation,  the  distribution  of  the  forced  ventilation- 
driven  flow  slabs  to  the  layers  of  the  receiving  room  and  their  rate  of 
extraction  from  the  layers  of  the  source  room  are  determined  with  use  of  the 
algorithm  FL0G02  discussed  in  Section  3.3.5. 


3.4.2.  Specification  of  the  Fan/Duct  System  Geometry  and  the  Geometry  and 
Fire  Environment  of  the  Connected  Room  Pairs 

Identify  the  two  rooms  of  the  facility,  or  the  room  and  the  outside 
environment,  connected  by  the  fan/duct  system  and  refer  to  the  Figure  3.4.1. 
Number  the  rooms  as  room  1 or  room  2 according  to  whether  it  is  connected  to 
the  fan/duct  system  inlet  or  outlet,  respectively.  Here,  the  system  inlet  and 
outlet  are  defined  as  the  side  of  the  fan/duct  flow  system  which  corresponds 
to  the  flow  inlet  and  outlet,  respectively,  under  normal  fan  operating 
conditions.  Normal  operating  conditions  lead  to  what  will  be  referred  to  here 
as  flow  in  the  normal  or  forward  direction,  from  room  1 to  room  2.  Flow  from 
outlet  to  inlet,  i.e.,  from  room  2 to  room  1,  will  be  referred  to  as  backward 
flow.  In  the  case  of  Option  1,  where  flow  direction  and  flow  rate  are 
specified,  the  number  identification  of  the  rooms  is  arbitrary. 

The  elevation  of  the  inlet  and  outlet  end-point  of  the  fan/duct  system  are 
characterized  by  the  elevation  of  their  midpoints,  ysysTEM  i>  room  i = 1 or 

2,  above  a datum  elevation.  The  vent  openings  at  the  duct  end-points  are 
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assumed  to  be  rectangular  with  sides  aligned  vertically  and  horizontally.  , 
the  depth  of  these  vents  in  room  i,  are  also  specified.  These  values  are 
required  in  order  to  determine  the  effect  on  the  fire  environment  of  the 
system  flow  at  instants  of  time  when  flow  through  the  system  is  drawn 
partially  from,  or  supplied  partially  to,  both  the  upper  and  lower  layer  of 
room  i . 


The  geometry  of  the  two  rooms  and  the  state  of  the  two -layer  environment  are 
specified  as  in  Section  3.3.2.  The  required  room  geometry  characteristics  are 
the  elevations  above  the  datum  elevation  of  the  floors  and  ceilings,  ypLooR  i 
and  Yceil  i’  ^respectively . 


The  instantaneous  states  of  the  room  environments  are  specified  by  the 
elevation  of  the  interfaces,  yL^YER  i’  separating  the  upper  and  lower  layers 
and  by  the  density,  absolute  temperature,  mass  concentration  mass  of  0^  , and 

mass  concentrations  of  the  other  (Npj^Qp  - 1)  products  of  combustion  in  the 
upper  and  lower  layers  of  room  i.  The  latter  are  designated  by  . [p^  ^1, 

'^U  , i ^0  2,U,i  ^0  2 , L , i ^ , U , i ^ , L , i ^ ^ ^ ^ ^PROD’ 

respectively,  for  the  upper  [lower]  layer.  At  instants  of  time  when  there  is 

no  upper  layer  [lower  layer],  i.e.,  when  yLAYER.i  = Yceil,!  [ytAYER.i  = 
ypLooR  i]>  upper-layer  [ lower- layer ] properties  are  not  used  in  the 
calculation.  Also,  as  in  Sections  3.2  and  3.3,  O2  and  other  product 
concentrations  would  be  required  only  when  the  simulation  under  consideration 
is  actually  being  used  to  predict  these  concentrations. 


The  final  specification  of  the  environment  is  SppLooR  i > difference 

between  the  absolute  hydrostatic  pressure  at  the  floor  of  each  room  and  some 
fixed,  absolute,  hydrostatic  datum  pressure,  Pdatum- 

Remarks  at  the  end  of  Section  3.3.2,  relative  to  situations  where  one  of  the 
"rooms"  connected  by  a natural  vent  is  the  outside  environment,  are  also 
applicable  here.  Also,  as  in  Sections  3.2  and  3.3,  all  above-mentioned 
elevations  are  relative  to  the  same  datum  elevation. 


3.4.3.  Specification  of  Fan/Duct  System  Characteristics 

As  mentioned  above,  there  are  two  options  for  specifying  the  fan/ duct  system 
characteristics.  In  Option  1,  the  flow  direction  and  volume  flow  rate  through 
the  system  are  specified.  By  electing  to  use  this  option,  one  assumes  that 
during  the  course  of  a simulated  fire  scenario  the  rate  of  volumetric  flow 
through  the  fan/duct  system  being  modeled  will  not  vary  significantly  from  a 
specified  design  value,  . Both  the  direction  and  the  magnitude  of  the 

volume  flow  rate  are  specified  by  the  value  of  . If  > 0 [Vpy^.j,  < 

0]  then  the  flow  is  forward  [backward]  flow  from  room  1 to  room  2 [room  2 to 
room  1 ] . 

When  using  Option  2,  the  direction  and  magnitude  of  the  flow  rate  through  the 
fan/duct  system  are  unknown.  To  carry  out  the  required  calculation  these  must 
be  determined  from  specified  fan/duct  characteristics  and  knowledge  of  the 
above-mentioned  parameters  which  define  the  two  room  environments. 


39 


I-Jhen  using  Option  2,  the  specified  fan/duct  characteristics  required  are  the 
fan  curve  and  the  duct's  flow  resistance. 

The  fan  curve  describes  the  operating  characteristics  of  the  fan.  It  is  a 
plot  of  the  cross-fan  pressure  difference,  App , as  a function  of  the  fan 
volume  flow  rate,  Vp  App  is  defined  as  the  difference  between  the  pressures 
at  the  fan's  outlet  and  the  inlet  sections.  This  is  identical  to  the 
difference  between  5pp  2>  5pp  , the  outlet  and  inlet  pressures, 

respectively,  as  measured  relative  to  the  datum  pressure. 


App  5pp  ^ 2 ■ ^Pf , 1 


(3.4.1) 


A generic  fan  curve  is  depicted  in  Figure  3.4.2. 

The  fan  curve  is  defined  analytically  by  linear  interpolation  between  and 
linear  extrapolation  beyond  two  or  more  pairs  of  specified  values,  [Vp  j, 

App  j] , J =1  to  JTERM  > 1,  which  identify  points  on  the  curve.  These  points 
indicate  the  volume  flow  rate  delivery  of  the  fan,  Vp  j,  at  discrete  cross-fan 
pressure  differences,  App  j. 

In  specifying  the  [Vp  j,  App  j]  data,  the  present  model  requires  the  App  j to 
be  monotonically  decreasing  with  increasing  Vp  j . Points  on  the  fan  curve  of 
Figure  3.4.2  would  satisfy  such  a requirement. 

The  flow  resistance  through  the  duct  would  be  dependent  generally  on  flow 
direction.  Designate  the  resistance  as  and  R2  for  forward  and  backward 
flow,  respectively.  Also  designate  the  density  of  the  flow  through  the 
fan/duct  system  as  will  be  determined  below,  this  will  also  be 

dependent  generally  on  flow  direction.  Then,  given  a non- zero  pressure  drop 
in  the  duct,  Ap^yj,^  , for  a mass  flow  rate  through  the  system  duct,  Mduct  ’ 
corresponding  volume  flow  rate,  the  non-zero  resistances,  R^  , are 

defined  as 


Ri  I APpjJ^  J j 1 / 2 ^ ^ ^ 2^1  ^ I ^^/(PQUQ^VQyj,.p) 

(3.4.2) 

4PduCT  ~ ('^PsYSTEM,!  * ^Pp  , 1 ^ (^Pf,2  ' ^PsYSTEM,2^ 

PdUCT  YSTEM  , 1 ‘ ySYSTEM,2) 


App  + (5psYSTEM . 1 


'^PSYSTEM,2)  PdUCT  S(ysYSTEM  , 1 ' ySYSTEM,2) 


where,  for  forward 
backward  flow,  i = 

'^PsYSTEM,2 

at  the  elevations. 


flow,  i = 1,  =1,  and  Appuf,^  would  be  > 0 , and  for 

2,  S2  = -1,  and  Apguf,^  would  be  < 0.  Also,  ^Psystem  1 
pressures  in  rooms  1 and  2,  relative  to  the  datum  pressure, 
ysYSTEM  1 ysYSTEM  2 ’ respectively,  of  the  inlet  and 
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outlet  system  endpoints.  Finally,  consistent  with  the  sketch  of  Figure  3.4.1, 
Eq . 3.4.2  assumes  that  there  is  a negligible  difference  in  elevation  between 
the  inlet  and  outlet  of  the  fan. 

If  the  duct  flow  resistance  is  negligible,  i.e.,  the  duct  resistances  can  be 
approximated  by  = ^2  = then  from  the  first  of  Eq . (3.4.2)  would 

also  be  approximated  as  zero  and  from  the  second  of  Eq . (3.4.2)  it  follows 
that 


if  = R2  = 0: 

^Pf  ~ '^^PsYSTEM.l  " ^PsYSTEM,2^  ' U C T S ( Ys Y S T E M , 1 ' YsYSTEM.Z^ 

(3.4.3) 


3.4.4.  First  Stage  of  the  Calculation:  Characteristics  of  the  Flow  Entering 
and  Leaving  the  Duct 

As  noted  in  Section  3.4.3,  for  Option  2 a priori  knowledge  of  the  direction  of 
the  flow  through  a particular  fan/duct  system  being  simulated  is  unavailable. 
It  can  be  from  room  1 to  room  2 or  from  room  2 to  room  1.  The  latter 
possibility  would  indicate  the  backward  flow  condition  where  flow  direction 
was  opposite  to  that  of  normal  fan  operation.  This  must  occur,  for  example, 
if  there  is  a fire  in  room  2,  the  outlet  room,  and  if  all  vents  and  other 
significant  penetrations  from  room  2 are  closed  except  for  the  room-2  end- 
point of  the  fan/duct  system. 

A priori  knowledge  of  flow  direction  depends  on  which  of  the  two  above- 
mentioned  algorithm  options  is  implemented.  To  keep  general  the  analysis  in 
the  discussion  to  follow,  the  subscripts  i used  in  the  variables  identified  in 
Section  3.4.2  are  replaced  by  iFR,  to  denote  the  source  room  from  which  the 
duct  flow  is  extracted,  and  by  iTO,  to  denote  the  receiving  room  into  which 
the  duct  flow  is  deposited.  Then,  one  task  of  the  first-stage  calculation  is 
to  determine  whether  iFR  = 1 and  iTO  = 2 (forward  flow,  or  normal  fan 
operation  with  > 0)  or  iFR  = 2 and  iTO  = 1 (backward  flow  through  the 

fan,  with  < 0)  . 


3. 4. 4.1.  Characterizing  the  Entering  Duct  Flow 
Consider  room  iFR. 

When  the  layer  interface  is  not  between  the  top  and  bottom  of  the  duct  opening 
a single  uniform  flow  slab  enters  the  duct  with  properties  identical  to  the 
specified  properties  of  the  lower  layer,  if  Yl^yer  ifr  — Ysystem  ifr  ^ifr/^’ 
the  upper  layer,  if  Ylayer.ifr  - Ysystem, ifr  ' ^ifr/^- 

When  the  layer  interface  is  between  the  top  and  bottom  of  the  duct  opening 
special  considerations  are  required. 
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Figure  3.4.3  is  an  illustration  of  the  flow  entering  the  duct  in  room  iFR  when 
the  layer  interface  is  between  the  top  and  bottom  of  the  vent  opening.  Define 
SpoucT  pressure  relative  to  the  datum  pressure  at  a location  somewhat 

inside  the  duct,  upstream  of  significant  duct- flow  development  and  mixing 
between  the  upper-  and  lower- layer  flow  contributions. 

Assume  that  the  flow  through  the  vent  and  into  the  entry  region  of  the  duct  is 
similar  to  the  flow  through  an  orifice  or  nozzle.  Assume  also  that  is 

small  enough  so  that  the  absolute  value  of  the  hydrostatic  pressure  difference 
over  the  height  of  the  duct  is  negligible  compared  to  the  absolute  value  of 
the  pressure  drop,  IAPentryI  ""  I^Pduct  ' '^Psystem  ifrI’  across  the  duct 
entrance  region.  Then,  from  Bernoulli's  equation,  e.g.,  Eq . (3.3.10)  for  e = 
0,  the  velocity  of  the  upper  and  lower  layer  flow-slab  components  of  the  flow 
in  the  entrance  region  of  the  duct  are  both  proportional  to  I Apg  | ^ ^ ^ . In 
particular,  the  mass  flow  rate  of  the  two  flow  slabs  from  the  upper  and  lower 
layer,  Mslabu.ifr  ^slabl.ifr-  respectively,  are  given  by 


■s  L A B U , 

, i FR  hCphy  ^ 

, i FR  ( 

1 ^Pentry 

Ipu, 

\l  ! 2 
iFR'' 

•S  L A B L , 

iFR  ~ hCg  hj^ 

, i FR  ^ 

1 ^Pentry 

IPl, 

Nl  / 2 
iFR'' 

In  the  above,  b is  the  width  of  the  vent  and  Cp  is  the  flow  coefficient. 

Also,  as  illustrated  in  Figure  3.4.3,  hy  ^ [h^  is  the  absolute  value  of  the 
difference  in  elevation  in  room  i between  the  top  [bottom]  of  the  vent  and  the 
layer  interface . 


. i “ ^SYSTEM,!  ■*"  ' YlAYER,! 

(3.4.5) 

^ , i ""  ytAYER.i  * ^ySYSTEM.i  ‘ /^ ) 

From  Eq.  (3.4.4)  it  follows  that 

^SLABU  , IFr/^SLABL  , iFR  ^ , i F R /^L  , i F R ^ ^ , i F R /^L  , i F R ^ ^ ^ (3.4.6) 


Assume  that  once  in  the  duct  the  two  streams  mix  adiabatically  to  form  a 
uniform-property  duct  flow  of  temperature  and  density  Pquct  where  the 

upper-  and  lower- layer  inlet  flows  can  both  be  modeled  as  the  same  perfect 
gas,  i.e. , 


PdUCt'^DUCT  Pu  , i Fr'^U  , i FR  , i F R , i F 1 


(3.4.7) 


Adiabatic  mixing  requires  that 
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M T 

“SLABU  , i FR  -^U  , i FR 


M T = I M IT 

^^SLABL  , iFR -^L  , iFR  I ‘‘■D  U C T I -^D  U C T 


(^SLABU.iFR  ^SLABL.iFR^'^DUCT 

(3.4.8) 


Now  Eqs . ( 3 . 4 . 6 ) - ( 3 . 4 . 8 ) lead  to 


^DUCT  ^DUCT , iFR 


■u  , i F R /^L  , 

, iFR ) , 

i Fr/^L  , 

')  1 / 2 1 / 

, i FR  ^ -1  / 

[1  + (hu, 

i F R /^L  , i 

FR ) , 

/d 

,iFR/^U,iFR^ 

(3.4.9) 


Thus,  if  the  flow  direction  (i.e.,  the  value  of  iFR)  is  known,  the  above  model 
for  flow  entry  into  the  duct  yields  an  estimate  of  Pduct’IFR  without  specific 
knowledge  of  the  value  of  or  its  entry  flow-slab  components,  I^slabu  ifr 

and  Mslabl  ifr-  Equation  (3.4.9)  will  be  used  when  the  present  forced-vent 
algorithm  is  implemented  in  the  option-2  mode  (see  Section  3.4.5).  At  such 
times  the  values  of  p^uct  '^rider  conditions  of  forward  and  backward  flow 
through  the  ventilation  system  (i.e.,  for  both  possibilities,  iFR  = 1 or  2) 
are  required  to  determine  the  direction  and  rate  of  flow  through  the  fan/duct 
system. 

Once  the  solutions  to  Pquct  determined  by  Eq.  (3.4.9),  Eq . (3.4.7)  leads 

'^DUCT 

'^DUCT""^U,iFR'^U,iFR/^DUCT  (3.4.10) 


Also,  once  the  solution  to  is  obtained  from  considerations  in  Section 

3. 4. 4. 3,  Eqs.  (3.4.6),  (3.4.8),  and  (3.4.10)  lead  to  ^ and  ^isLABL.iFR 


K 


SLABL.iFR 


UC  T 1 UC  T / ^ '^U  , i FR  [ , i Fr/^L  , i FR  ^ , i Fr/^L  , i FR  ) ^ ^ 

, i FR  ^ 

(3.4.11) 


M, 


SLABU , iFR 


INDUCT  I ‘ ^SLABL.iFR 


1 / 2 


INDUCT  I'^DUCt/^'^L,  IFr[^^L,  IFr/^U  , IFR^^^L,  IFr/  ^U,  IFR^ 

+ T 


U , I FR 

(3.4.12) 
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3. 4. 4. 2.  Determining  the  Entering  Duct  Flow  When  the  Flow  Direction  and  Flow 
Rate  Are  Specified 


When  using  Option  1,  the  values  of  and  iFR  are  known,  and  the  values  of 

Pduct  ’ "^DUCT  ’ ^DUCT  entering  components  of  the  duct  flow  slabs  are 

determined  according  to: 


If  the  interface  is  between  the  top  and  bottom  of  the  room- iFR  vent: 

There  are  two  flow  slabs  entering  the  vent;  Pquqj  and  are 

determined  from  Eqs . (3.4.9)  and  (3.4.10),  respectively; 

Mj^uct  ^ ^DUCT^DUCT’  , using  the  latter  value  for  , 

Mslabl  iFR  ^SLABU  iFR  determined  from  Eqs.  (3.4.11)  and 

(3.4.12),  respectively. 

(3.4.13) 

If  the  interface  is  at  or  below  the  bottom  of  the  room- iFR  vent: 

vent  from  the  upper  layer; 

^DUCT^DUCT’  ^SLABU,iFR  ^ 

(3.4.14) 


There  is  one  slab  of  flow  entering  the 

T. 


^QUCT  , iFR  > -^DUCT 


IM, 


DUCT 


and  M. 


= T • M 

■^u  , iFR  ’ “duct 


SLABL.iFR 


= 0. 


If  the  interface  is  at  or  above  the  top  of  the  room- iFR  vent: 


There  is  one  flow  slab  entering  the  vent  from  the  lower  layer; 

PqUCT  ^ ^L.iFR’  "^DUCT  ^ "^L.iFR’  ^DUCT  ^ ^DUCT^DUCT’  ^SLABL.iFR 

IM. 


^DUCT 


and  Mg^ABu.iFR  ^ 


(3.4.15) 


From  the  above  and  from  other  specified  information,  all  other  properties  of 
the  entering  flow  slabs  and  the  duct  flow  can  also  be  determined. 


3. 4. 4. 3.  Determining  the  Entering  Duct  Flow  When  the  Flow  Rate  and  Flow 
Direction  Are  Unknown 

In  the  Option- 2 mode,  it  is  necessary  to  determine  the  entering  flow  slab 
components  of  the  duct  flow  and  the  properties  of  the  duct  flow  itself  as  in 
Option  1.  Here,  the  system  flow  rate  and  the  flow  direction  (i.e.,  the  value 
of  iFR)  are  not  known. 

To  obtain  the  flow  rate  it  is  necessary  first  to  determine 
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'^PSYSTEM 


= 5p, 


SYSTEM  , 1 


SYSTEM , 2 


(3.4.16) 


- 6p 


where  ^ and  ^Pg^g^E^  2 obtained  from: 


For  i = 1 and  2 : 


ySYSTEM.i  ^ yLAYER.i’ 

^PsYSTEM.i  ~ ^PfLOOR,! 


Pl , i S^^LAYER , i ~ ypLOOR.i^ 

Pu , i S^ys YSTEM , i ‘ yLAYER.i^ 

(3.4.17) 


ySYSTEM.i  - ytAYER.i- 

'^PsYSTEM.i  ~ ^PpLOOR,!  * , i § ^^ys  Y S T E M , i ' ypLOOR.i^ 

(3.4.18) 


Next,  the  two  possible  P^uct  values,  corresponding  to  forward  and  backward 
flow,  are  determined.  These  are  designated  as  P^uct  1 there  is  forward 
flow  (i.e.,  if  iFR  = 1)  and  Pquct  2 there  is  backward  flow  (i.e.,  if  iFR  = 
2) . The  values  of  P^uct  i determined  one  at  a time  according  to  the 

Option-1  determination  of  Pquct  pi^esented  in  Eqs  . ( 3 . 4 . 13 ) - ( 3 . 4 . 14)  . 

Finally  it  is  necessary  to  determine 


^ySYSTEM  ySYSTEM.l  ' ySYSTEM,2 


(3.4.19) 


Now,  with  the  values  of  Apj^gj^j.,  , 4y2Yg.pEM>  Pduct  i’  P^® 

characteristics  of  the  fan-duct  system,  namely,  the  R^ , JTERM  and  the  [Vp  j , 
App  j ] , J = 1 to  JTERM,  the  algorithm  described  below  (in  Section  3.4.5)  is 
used  to  determine  M^ucx • 

As  will  be  mentioned  in  Section  3.4.5,  there  are  situations  where  the  adopted 
model  of  the  fan-duct  system  can  yield  two  solutions  for  the  flow,  a forward 
flow  solution  and  a backward  flow  solution.  Under  such  a circumstance,  the 
solution  adopted  here  is  the  one  that  corresponds  to  the  flow  direction 
determined  at  the  previous  time  step  in  the  fire  simulation. 
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3. 4. 4. 4.  Characterizing  the  Flow  as  it  Leaves  the  Fan-Duct  System  and 
Penetrates  the  Receiving  Room 

The  first-stage  calculation  concludes  with  a characterization  of  the  known, 
uniform-property  duct  flow  as  it  leaves  the  fan/duct  system  and  penetrates  the 
layers  of  the  receiving  room,  room  iTO . In  terms  of  a flow  slab  description 
of  the  penetrating  flow,  the  following  determination  is  made; 


If  the  interface  in  room  iTO  is  between  the  top  and  bottom  of  the  room- 
iTO  vent : 

There  are  two,  identical -property  flow  slabs  leaving  the  vent. 

The  upper  flow  slab  consists  of  the  portion  of  the  duct  flow 
between  the  elevation  of  the  layer  interface  and  the  top  of  the 
vent,  and  the  lower  flow  slab  consists  of  that  portion  of  the  duct 
flow  between  the  interface  and  the  bottom  of  the  vent.  The  upper 
and  lower  flow  slabs  have  depths  of  hy  and  h^  ^.j.q  which  are 

determined  by  Eq . (3.4.5),  and  they  penetrate  the  upper  and  lower 
layers  of  room  iTO,  respectively.  The  mass  flow  rate  in  the  upper 
and  lower  slabs  are  designated  as  Mslabu  ifr  ^slabl  ifr’ 

respectively.  These  are  determined  from 

^SLABU.iFR  , i T 0 T 0 ^ I U C T I 

^SLABL.iFR  ^ , i T 0 T 0 ^ I U C T I I^DUCtI  " ^SLABU.iFR 

(3.4.20) 


If  the  interface  in  room  iTO  is  at  or  below  the  bottom  of  the  room- iTO 
vent ; 


There  is  one  flow  slab  leaving  the  vent  and  penetrating  the  upper 
layer . 


(3.4.21) 


If  the  interface  in  room  iTO  is  at  or  above  the  top  of  the  room- iTO 
vent : 


There  is  one  flow  slab  leaving  the  vent  and  penetrating  the  lower 
layer . 


(3.4.22) 


The  properties  of  the  entering  flow  slabs  are  those  of  the  duct  flow. 
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3.4.5.  The  Direction  and  Flow  Rate  Through  a Fan/Duct  System 

This  Section  describes  the  basis  for  calculating  the  flow  rate  through  a 
fan/duct  ventilation  system  of  the  type  depicted  in  Figure  3.4.1.  In  carrying 
out  the  calculation,  modeling  of  the  system  follows  an  extension  to  the  ideas 
presented  in  reference  [3.4.1]. 

In  Figure  3.4.1,  Ysystem  i’  ysYSXEM  2 mid- elevation , relative  to  the 

datum  elevation,  of  the  system  end-points  in  rooms  1 and  2,  respectively. 

Also,  5p2Y2TEM  1’  ^PsYSTEM  2 pressures  at  these  end-points,  relative 

to  the  datum  pressure.  As  defined  in  Section  3.4.2  and  as  indicated  in  the 
figure,  the  end-points  1 and  2 correspond  to  the  inlet  and  outlet  of  the 
system  under  normal  fan  operating  conditions. 

To  carry  out  the  present  calculation,  several  parameters  of  the  system  must  be 
specified.  These  are:  the  pressure  and  elevation  differences,  ApgysxEM 
AysYSTEM’  defined  in  Eqs . (3.4.16)  and  (3.4.19);  the  forward  and  backward 
flow  resistance  of  the  duct,  and  R2  , respectively,  as  defined  in  Eq . 

(3.4.2);  and  the  two  possible  densities  of  the  duct  flow,  designated  as 
Pduct  1 ’ solution  flow  is  a forward  flow  from  the  inlet  room  1 to  the 

outlet  room  2,  and  Pquct  2>  solution  flow  is  a backward  flow  from  the 

outlet  room  2 to  the  inlet  room  1. 

Also  required  for  the  calculation  are  two  or  more  pairs  of  specified  values 
[Vp  j,  App  j ] , J = 1 to  JTERM  > 1,  which  identify  points  on  the  fan-operating 
curve.  As  defined  in  Section  3.4.3,  the  Vp  j's  are  the  volume  flow  rates 
delivered  by  the  fan  under  the  corresponding  cross -fan  pressure  differences 
App  J,  where  App  is  defined  in  Eq . (3.4.1).  As  indicated  in  Section  3.4.3,  in 
simulating  fan  operation  the  fan  curve  is  approximated  by  linear  interpolation 
between  and  linear  extrapolation  beyond  the  specified  JTERM  fan-curve  data 
points.  As  indicated  in  Figure  3.4.2,  the  fan  curve  must  be  such  that  App  is 
monotonically  decreasing  with  increasing  Vp . In  terms  of  the  fan  curve, 
forward  flow  is  associated  with  positive  values  of  the  volume  flow  rate,  Vp , 
while  backward  flow  is  associated  with  negative  values  of  Vp . 

Let  the  simulated  fan  curve  be  defined  as 


Vp  = f(App) 


(3.4.23) 


and  define  its  Vp  = 0 intercept  as  App  q,  i.e., 

f(App  0)  = 0 (3.4.24) 

Also,  using  the  above  - specif ied  parameters  and  referring  to  Eqs.  (3.4.2), 
(3.4.16),  and  (3.4.19),  define  and  evaluate  the  adjusted  pressure  differences, 

^Padj, i - 
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SYSTEM 


^DUCT , i S^^S YSTEM 


(3.4.25) 


^Pa  D J , i 


Ap 


where  i = 1 or  2 for  forward  or  backward  flow,  respectively. 

Using  Eq . (3.4.25)  in  Eq . (3.4.2)  leads  to  the  following  result  for  the  duct 
mass  and  volume  flow  rates: 


^DUCT.l  Pd  UCT,  INDUCT,!  I ^Pp  ^Pa  D J , 1 I ^ ^ ^ (^Pp  '*■  ^PaDJ.I^  ^ ® 

(3.4.26) 

^DUCT,2  ^ PdUCT,2^DUCT,2  ^ " I ^Pp  ^Pa  D J , 2 I ^ ^ ^ /^2  (^Pp  "*■  ^PaDJ,2)  ^ 

(3.4.27) 


For  a valid  fan/duct  system  flow  solution,  the  duct  volume  flow  rate  in  Eqs . 
(3.4.26)  or  (3.4.27)  and  its  corresponding  mass  flow  rate  must  be  identical  to 
the  fan  volume  flow  rate  of  Eq . (3.4.23)  and  its  corresponding  mass  flow  rate. 
Equating  these  volume  flow  rates  leads  to  the  following  result: 


For  non- zero  R. , App  solutions  to 


|App  + Ap^p  j I ^ ^ ) 

or 

I ^Pp  ^Pa  DJ,2  I^^^/*'PdUCT,2^2^ 


f (App ) if  f (App  ) > 0 

(3.4.28) 

f (App ) if  f (App ) < 0 


and  the  corresponding  = Vp  solutions  that  would  be  obtained  from 

Eq . (3.4.23),  would  represent  valid  [Vp , App]  solution  pairs  to  the 
desired  fan/duct  flow  problem. 


If  the  duct  flow  resistance  is  negligible  then  Eqs.  (3.4.3)  and  (3.4.25)  lead 
to  the  following  alternative  to  the  result  of  Eq . (3.4.28): 
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For  = ^2  ^ ^Pf  solutions 


^Pf  = -^Padj.i  if  f(-^PADJ,i)  > 0 

or  (3.4.29) 

App  = 2 if  f(*^PADJ,2^  ^ i* 

and  the  corresponding  ^ ^f  solutions  that  would  be  obtained  from 

Eq . (3.4.23)  would  represent  valid  [Vp , App ] solution  pairs  to  the 
desired  fan/duct  flow  problem. 


If  Eqs . (3.4,28)  or  (3.4.29)  do  not  lead  to  any  solutions,  then  the  solution 
to  the  duct/fan  flow  problem  is  one  of  zero  flow,  i.e.. 


'DUCT 


V. 


0 and  App  = Ap 


F , 0 


if  there  are  no  solutions  to 
Eqs.  (3.4.28)  or  (3.4.29). 


(3.4.30) 


Figure  3.4.4  is  a graphical  presentation  of  the  problem  of  solving  Eqs. 
(3.4.28)  or  (3.4.29).  Plotted  there  is  the  universal  fan  curve  of  Eq . 
(3.4.23),  the  solid  line,  and  in  the  dashed  lines,  parametric  plots  of  Eq . 
(3.4.26)  for  forward  duct  flow  on  the  right  side  of  the  Ap  axis,  and 
parametric  plots  of  Eq . (3.4.27)  for  backward  duct  flow  on  the  left  side  of 
the  Ap  axis. 

The  parametric  curves  on  the  right  side  of  Figure  3.3.4  characterize  forward 
flow  through  the  duct.  According  to  Eq . (3.4.26),  these  correspond  to  plots 


of  App  vs  V for  different  values  of  Ap 


AD  J , 1 


r(F  ) 


r(F) 


r(F) 


0, 


r(  F ) 


etc . 


At  any  instant  and  for  any  duct  design,  the  value  of  Ap^^j  ^ calculated  from 
Eq . (3.4.25)  would  be  equal  to  one  of  these  values  (or  to  an  interpolated 

value) . Then  the  plot  of  App  vs  V for  the  duct  would  be  identical  to  the 
corresponding  curve  shown  in  the  figure  (or  to  an  interpolated  curve 

between  the  ones  shown) . 


In  an  analogous  way,  the  parametric  curves  on  the  left  side  of  Figure  3.3.4 
characterize  backward  flow  through  the  duct.  According  to  Eq.  (3.4.27),  these 
correspond  to  plots  of  App  vs  V for  different  values  of  Ap^^j  2 = 


r(B  ) 


r(  B ) 


0,  etc.  At  any  instant  and  for  any  duct  design,  the  value  of 


calculated  from  Eq . (3.4.25)  would  be  equal  to  one  of  these 


cr(B) 

^Pad  J , 2 

values  (or  to  an  interpolated  value) . Then  the  plot  of  App  vs  V for  the  duct 
would  be  identical  to  the  corresponding  curve  shown  in  the  figure  (or  to 

an  interpolated  curve  between  the  ones  shown) . 

Note  that  if  duct  flow  resistance  is  zero,  then  from  Eq . (3.4.29)  the 
parametric  duct- flow  curves  become  the  parametric  horizontal  lines 
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> 0 and  App  = -APadj , i 


0 


(3.4.31) 


V 


DUCT 


if  = 


for  forward  flow,  and 


^DUCT  ^ and  App  ■'^Padj,2  ^^^2  ^ (3.4.32) 

for  backward  flow. 

In  Figure  3.4.4,  a valid  forward- flow  solution  would  correspond  to  the  unique 
intersection  of  the  solid  fan  curve  and  the  dashed  duct  curve  on  the  right 
which  corresponds  to  the  particular  value  of  = Ap^pj  ^ calculated  in  Eq . 

(3.4.25).  In  the  figure  it  can  be  seen  from  the  intercepts  of  the  ^ curves 
and  the  App  axis  that  a forward-flow  solution,  App  soln > exists  if  and  only  if 
^PpDJ.i  satisfies  ^ = '^Padj  x ^Pf  , 0 • 

Analogous  to  the  forward  flow  solution,  a valid  backward- flow  solution  would 
correspond  to  the  unique  intersection  of  the  solid  fan  curve  and  the  dashed 
duct  curve  on  the  left  which  corresponds  to  the  particular  value  of  = 

APadj  2 calculated  in  Eq . (3.4.25).  In  the  figure  it  can  be  seen  from  the 
intercepts  of  the  curves  and  the  App  axis  that  a backward- flow  solution, 

App  soln > exists  if  and  only  if  Ap^pj  2 satisfies  -a^®^  = -Apppj  2 > ^Pp  o- 

The  graphical  presentation  of  Figure  3.4.4  clarifies  the  fact  that  there  are 
four  possible  classes  of  solution  to  Eqs . (3.4.28)  or  (3.4.29),  namely:  a 
single  forward-flow  solution  and  no  backward- flow  solution,  e.g.,  Ap^pj  ^ = 
and  Ap^pj  2 = ^3®^  [see  Figure  (3. 4. 5. a)],  a single  backward- flow 
solution  and  no  forward-flow  solution,  e.g.,  Ap^pj  ^ ^Padj  2 ^ 

ct[®^  [see  Figure  (3.4.5.b)],  one  forward-flow  solution  and  one  backward- flow 
solution,  e.g.,  Ap^pj  ^ = *^2^^  '^Padj  2 ^ [see  Figure  (3.4.5.c)],  and 

no  solutions  at  all,  e.g.,  Ap^pj  ^ ^ ^Padj  2 *^2^^  [see  Figure 

(3.4.5.d)].  In  the  first  two  classes,  the  unique  solution  obtained  is  taken 
to  be  the  desired  solution  to  the  fan/duct  flow  problem.  In  the  latter  class, 
the  solution  to  the  fan/duct  flow  problem  is  taken  to  be  that  of  Eq . (3.4.30). 
As  mentioned  at  the  end  of  Section  3. 4. 4. 3,  when  both  forward-  and  backward- 
flow  solutions  are  obtained  (i.e.,  the  third  class  of  solution)  the  solution 
to  the  fan/duct  flow  problem  is  taken  to  be  the  solution  with  a flow  direction 
identical  to  the  flow  direction  determined  at  the  previous  instant  of  time  in 
the  fire  simulation. 

The  actual  solution  to  Eq . (3.4.28)  makes  use  of  the  root- finding  algorithm 
RTSAFE  and  its  associated  subroutine  presented  in  [3.4.2].  The  RTSAFE 
subroutine  is  a part  of  the  FANRES  subroutine  of  CCFM. VENTS  which  is  included 
in  Part  III  [1.1.2]  of  this  work. 
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3.4.6.  Second  Stage  of  the  Calculation:  Distribution  of  the  Duct  Flow  Slabs 
to  the  Receiving  Layers  and  Rate  of  Extraction  From  the  Source  Layers 

The  result  of  the  first-stage  calculation  of  Section  3.4.4  is  a description  of 
the  flows  in  and  out  of  the  fan/duct  system  where  these  flows  are  described  in 
the  manner  of  the  VENTHP  natural  vent  flow  slabs  of  Section  3.3.4.  In  the 
second-stage  calculation,  the  distribution  of  the  present  forced  ventilation- 
driven  flow  slabs  to  the  layers  of  the  receiving  room  and  their  rate  of 
extraction  from  the  layers  of  the  source  room  are  determined  according  to  the 
algorithm  FL0G02  discussed  in  Section  3.3.5. 

As  will  be  recalled,  FL0G02  provides  the  contributions  to  the  rates  of  change 
of  the  upper  and  lower  layers  of  two  immediately  adjacent  rooms  as  a result  of 
the  one-to-six  slabs  of  uniform-property  flow  which  pass  through  a vertical 
vent  in  their  common  wall  segment. 

In  the  configuration  of  Figure  3.4.1,  the  flow  slabs  associated  with  the  vent 
flows  at  either  end  of  the  fan/duct  system  have  been  determined;  the  flow  from 
the  source  room  to  the  duct,  at  one  end,  and  the  flow  from  the  duct  to  the 
receiving  room,  at  the  other.  Two  successive  applications  of  the  FL0G02 
algorithm  are  used.  In  the  first  application,  the  pair  of  immediately 
adjacent  rooms  identified  in  the  algorithm  is  the  source  room,  room  iFR,  and  a 
"dummy  room,"  associated  with  the  outside  environment.  The  latter  is  taken  as 
a surrogate  for  the  FL0G02  receiving  room.  In  the  second  application,  the 
pair  of  adjacent  rooms  is  the  receiving  room,  room  iTO , and  a "dummy  room," 
now  taken  as  a surrogate  to  the  FL0G02  source  room.  As  mentioned  at  the  end 
of  Section  3. 3. 5. 6,  when  one  of  the  rooms  of  a room  pair  in  a FL0G02 
calculation  is  a space  associated  with  the  outside  environment,  modifications 
to  that  outside  space  as  a result  of  leaving  or  entering  flow  slabs  are 
ignored.  Therefore,  in  determining  the  action  of  the  forced  ventilation 
system  the  result  of  the  two  applications  of  FL0G02  are  1)  to  calculate  the 
rates  of  loss  of  mass,  enthalpy,  and  products  from  the  layers  of  the  source 
room;  and  2)  to  calculate  the  rates  of  addition  of  mass,  enthalpy  and  products 
to  the  layers  of  the  receiving  room. 

The  contributions  of  forced  ventilation  flows  to  the  layers  of  all  rooms  of  a 
facility  is  obtained  by  carrying  out  the  above  fan/duct  calculations 
sequentially  for  all  fan/duct  systems  of  the  facility.  For  each  system  this 
would  involve  both  stages  of  the  calculations  of  Section  3.4. 


3.4.7.  The  Algorithms  and  Associated  Subroutines,  VENTF,  FL0G02,  and  FANRES , 
Used  to  Carry  Out  the  Two  Stages  of  the  Forced  Vent  Flow  Calculation 

The  methodology  described  in  Sections  3.4.1  - 3.4.6  for  estimating  the  rates 
of  exchange  of  mass,  enthalpy,  and  products  of  combustion  through  forced- 
ventilation  fan/duct  systems  is  implemented  in  three  separate  algorithms  named 
VENTF,  FANRES,  and  FL0G02 . 

VENTF  implements  the  calculations  of  Stage  1 which  correspond  to  the 
presentation  in  Section  3.4.4.  FL0G02  implements  the  calculations  of  Stage  2 
which  correspond  to  the  presentation  in  Sections  3.4.6  and  3.3.5.  In  cases 
where  the  flow  through  the  fan/duct  system  is  determined  from  the  duct  flow 
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resistance  and  the  fan  operating  curve,  the  algorithm  FANRES  is  used  in  the 
VENTF  calculation.  The  calculation  of  FANRES  corresponds  to  the  presentation 
in  Section  3.4.5. 

The  use  of  the  three  algorithms  of  the  forced  vent  calculation  is  depicted  in 
Figure  3.4.6. 

A concise  presentation  of  VENTF,  FANRES,  and  FL0G02  and  their  associated 
FORTRAN  subroutines  appears  in  Part  111  [1.1.2]  of  this  work. 
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4.  Specifying  the  Problem  for  a Particular  CCFM. VENTS  Fire  Simulation 

To  carry  out  a particular  fire  simulation  with  CCFM. VENTS  the  user  must  supply 
information  on  the  characteristics  of  the  rooms,  the  natural  vents,  the  forced 
vents,  the  fires,  and  the  outside  environment.  Also  required  are  the  initial 
conditions  which  describe  the  state  of  the  two -layer  environment  in  each  room 
of  the  simulated  facility.  The  required  inputs  will  be  discussed  here  briefly 
Other  required  inputs,  not  associated  with  the  physical  aspects  of  the 
simulation,  will  be  discussed  in  Part  II  [1.1.1]  of  this  work. 

The  detailed  mechanics  of  implementing  all  CCFM. VENTS  input  are  described  in 
user's  reference  guide  of  Part  IV  [1.1.3]. 


4.1.  Room  Geometry,  Vent  Characteristics,  and  the  Outside  Environment 

Independent  sets  of  positive  integers  are  used  to  number  and  designate  the 
inside  rooms,  the  natural  vents,  and  the  forced  vents  (i.e.,  the  fan/vent 
systems)  of  the  facility. 

As  will  be  discussed  below,  when  dealing  with  wind  environments  and  for  other 
applications,  it  will  be  useful  to  divide  conceptually  the  outside  environment 
into  a number  of  different  spaces.  These  outside  spaces  would  also  be  given 
designations  from  an  independent  set  of  positive  numbers.  When  identifying 
one  of  these  outside  spaces,  e.g.,  to  specify  that  it  is  one  of  the  two  spaces 
joined  by  a vent,  the  assigned  integer  for  the  outside  space  is  distinguished 
from  that  of  the  inside  room  by  use  of  the  negative  of  its  designated  integer 
value.  In  what  is  likely  to  be  the  most  typical  of  fire  scenarios,  when  the 
outside  environment  is  taken  to  be  quiescent  and  uniform,  only  one  outside 
space  is  required.  This  space  is  designated  as  outside -environment  space  1. 

The  specified  geometric  properties  of  each  room  i are  ypLooR  i height 

and  area. 

The  specified  properties  of  each  natural  vent  i are:  the  numerical  designation 
of  the  two  Inside  rooms  (or  the  Inside  room  and  outside  space)  connected  by 
the  vent,  one  of  these  being  identified  arbitrarily  as  the  "from"  room  and  the 
other  as  the  "to"  room;  the  area  of  the  vent;  and  the  elevations,  above  the 
floor  of  the  "from"  room,  of  the  bottom  and  top  of  the  vent. 

The  specified  properties  of  each  forced  vent  i are:  the  numerical  designation 
of  the  "from"  room  and  the  "to"  room  under  normal  fan-operating  conditions; 
the  elevations  above  the  floors  of  their  respective  rooms  of  the  midpoints  of 
the  fan/duct  system  endpoints;  the  heights  of  the  rectangular  duct  end-point 
vents  in  the  "from"  and  in  the  "to"  rooms;  the  flow  resistances  of  the  duct 
during  forward  and  backward  flow  conditions;  and  a set  of  pairs  of  values  of 
volume  flow  rate  and  cross -fan  pressure  which  identify  points  on  the  fan 
curve . 

The  specified  properties  of  each  outside  space  i are:  a reference  elevation, 
relative  to  the  datum  elevation,  where  ambient  properties  of  the  space  are 
specified;  the  reference  pressure,  relative  to  the  datum  pressure,  at  the 
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latter  elevation;  the  density  of  the  space;  and  the  concentration  in  the  space 
of  all  simulated  products  of  combustion. 


4.2.  The  Fire 

The  specified  properties  of  the  fires  are:  the  numbers  of  the  rooms  that 
contains  a fire;  the  elevation  of  the  fire  above  the  floor;  and  sets  of  pairs 
of  values  of  energy  release-rate  and  time  which  identify  points  on  the 
specified  curves  of  fire-energy-release-rates  vs  time.  In  the  simulation,  for 
each  of  the  specified  fires  the  fire-energy-release-rate  at  an  arbitrary  time 
is  obtained  by  interpolating  linearly  between  or  extrapolating  linearly  beyond 
the  specified  data  points. 


4.3.  The  Products  of  Combustion 

An  independent  set  of  positive  integers  is  used  to  number  and  designate  the 
different  products  of  combustion  whose  concentration  is  to  be  simulated.  At 
an  arbitrary  instant  of  time  it  is  assumed  that  the  fire's  generation  rate  of 
product  of  combustion  k,  Pj^  source  ’ proportional  to  the  simulated  energy- 
release-rate  of  the  fire,  i.e., 

^k, SOURCE  ^ factorj^ -Q  (4.3.1) 


Thus,  for  product  k,  the  specified  property  of  the  product  generation  rate  is 
the  constant,  factorj^  . The  default  product  of  combustion  1 is  oxygen  and, 
according  to  Eqs.  (3.2.2),  (3.2.3),  and  (4.3.1),  a default  value  for  factor^^ 
is  specified  as  -0 . 076 ( 10" ® ) (kg  of  0^) /^ ■ 


4.4.  Initial  Conditions 

The  environment  in  each  room  i of  the  facility  must  be  initialized.  The 
initial  properties  specified  are:  ^Pp^ooR  i’  elevation  of  the  layer  interface 
above  the  floor,  ^ , p^  ^ , and  u i L i each  product  of 

combustion,  k.  Note  that  if  the  upper  layer  is  initially  of  zero  thickness, 
then  py  . and  the  Pj^  ^ ^ can  be  specified  arbitrarily.  Similarly,  if  the 
upper  layer  initially  fills  the  room,  p^  ^ and  the  Pj^  l i specified 

arbitrarily . 

As  will  be  seen  in  Part  IV  [1.1.3],  a simple  procedure,  appropriate  for  use  in 
most  simulations,  is  available  for  initializing  all  of  the  above  variables. 

In  this  procedure:  1)  all  rooms  of  the  facility  are  initialized  with  a zero- 
thickness upper-layer;  and  2)  the  density  and  concentration  of  each  product  of 
combustion  are  taken  to  be  uniform,  at  specified  values,  throughout  all  rooms 
of  the  facility  and  throughout  all  outside  spaces.  With  the  above  specified 
information,  the  exact  values  of  ^Pp^ooR  i rooms  is  calculated 

automatically.  The  procedure  leads  to  the  desired  condition  of  initial-zero- 
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flow  between  all  inside  - to  - inside  and  inside- to-outside  natural  vents  of  the 
facility . 

Having  used  the  above  initialization  procedure,  the  user,  as  required,  can 
then  make  initialization  modifications  to  the  outside  spaces  and  to  selected 
inside  rooms. 


4.5.  Specifying  a Simulation  Which  Accounts  for  Conditions  of  Wind  Loading 

Except  in  unusual  situations  where  it  is  reasonable  to  model  a facility  as 
being  completely  sealed  off  from  the  outside  environment,  outside  wind 
conditions  can  play  a significant  role  in  determining  the  fire  - generated 
environment  within  the  facility.  This  is  because:  1)  the  wind  dynamics  affect 
significantly  the  exterior  pressure  distribution  local  to  the  outside- 
environment  -exposed  side  of  any  particular  natural  or  forced  vent;  and  2)  such 
an  exterior  pressure  distribution  plays  a major  role  in  determining  the 
characteristics  of  the  flow  through  the  vent,  e.g.,  this  distribution  would  be 
computed  and  used  in  the  algorithm  VENTHP . 

The  free-field  wind  velocity  relative  to  the  orientation  of  the  facility  and 
the  specific  location  on  the  exterior  of  the  facility  structure  will  determine 
the  local  pressure  distribution  there.  The  location  and  shape  of  nearby 
structures  can  also  play  a role  in  determining  the  exterior  pressure 
distribution. 

The  CCFM. VENTS  capability  of  designating  multiple  outside  spaces  can  be  used 
to  simulate  the  pressure  distributions  at  different  vent  locations  on  the 
exterior  of  a facility  in  a wind  environment.  One  possible  method  of  using 
CCFM. VENTS  to  simulate  such  wind  conditions  is  outlined  briefly: 

Specify  , the  density  of  a quiescent  environment,  5pj^j.p  , its  reference 

pressure  at  its  reference  elevation,  and  y^gp , its  reference  elevation. 
Designate  this  environment  as  outside  space  1. 

Now  assume  that  the  environment  is  not  quiescent,  but  that  it  has  a uniform 
far-field  wind  velocity  with  wind  speed, 

Consider  a vent  of  a facility  which  is  exposed  to  an  outside  wind  environment. 
Assume  that  an  analysis  of  wind  flow  around  the  facility  allows  one  to 
determine  the  local  pressure  coefficient,  Cpj^pgguRE'  location  of  the 

exposed-side  of  the  vent  of  interest  (see,  e.g.,  [4.5.1]).  Thus,  if 
^Pquiescent  ^PwiND  relative  outside  pressures  local  to  the  vent  in 

the  quiescent  and  wind  environments,  respectively,  then,  by  definition  of 
r 

'^PRESSURE  ’ 

'^PwiNO  ""  ^PqUIESCENT  *^PRESSURE  (^AMB^WIND^/2^  (^-5.1) 

In  view  of  Eq . (4.5.1),  a new  outside  space  is  introduced,  where  the  new  space 
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is  defined  as  having  the  same  properties  as  the  original  quiescent -environment 
space,  but  with  a new  reference  pressure  specified  to  be 


(Sp 


REFNEW 


= 5p 


REF 


*^PRESSURE  (^AMb’^WIND 


72) 


(4.5.2) 


The  "room"  on  the  exposed  side  of  the  vent  of  interest  is  now  designated  as 
this  new  outside  space. 

Now  the  above  procedure  is  repeated  for  each  exposed  vent  of  the  facility.  At 
the  end  of  the  process  one  outside  space  will  be  defined  for  each  unique  value 

*^PRESSURE  • 

With  the  modified  input  data,  the  desired  effect  of  wind  loading  would  be 
taken  account  of  in  the  resulting  CCFM. VENTS  simulation. 


4.6.  Stack  Effect 

Flows  through  natural  vents  located  on  the  exterior  walls  of  the  facility 
(e.g.,  cracks,  open  windows,  or  doors)  can  be  driven  by  cross-vent  pressure 
differences  which  are  set  up  as  a result  of  a density  difference  between  the 
interior  of  a facility  and  the  outside  environment.  The  flow  exchange 
phenomenon  is  called  stack  effect  [4.6.1]. 

The  phenomenon  of  stack  effect  can  be  very  important  in  relatively  tall 
buildings  during  winter  or  summer  conditions.  The  winter  condition  involves  a 
relatively-uniform  interior-building  environment  which  has  a lower  density 
(higher  temperature)  than  the  density  (temperature)  of  the  outside  ambient. 

The  summer  stack-effect  phenomenon  involves  a relatively  high  interior  density 
(low  interior  temperature)  and  relatively  low  outside  ambient  density  (high 
outside  temperature) . 

The  combined  phenomena  of  stack- effect- driven  flows  and  fire -driven  flows  can 
be  simulated  with  CCFM. VENTS.  One  possible  way  to  do  this  is  outlined 
briefly: 

Initialize  the  simulation  according  to  Section  4.4,  setting  the  uniform  inside 
and  outside  density  at  identical  values  corresponding  to  the  desired  inside 
density.  Now  modify  the  outside  environment  by  specifying  the  desired  outside 
density . 

With  the  modified  input  data,  the  CCFM. VENTS  simulation  will  include  the 
desired  phenomenon  of  stack  effect. 

Note  that  the  ideas  on  wind- load  simulation  of  Section  4.5  can  be  combined 
with  the  present  ideas  on  stack-effect  simulation.  Thi.s  leads  to  a dapability 
of  simulating  the  combined  effect  of  fire,  wind,  and  stack  effect. 
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5.  Summary  and  Overview  of  Parts  II -IV 


This  is  the  first  of  a four-part  report  which  documents  all  aspects  of 
CCFM. VENTS,  the  first-stage  prototype  application  of  the  CCFM  concept. 
Introductory  remarks  of  Section  1 discussed  the  generic  features  of  the  CCFM, 
and  the  specific  features  of  CCFM. VENTS. 

The  main  objective  of  this  Part  1 document  was  to  present  a comprehensive 
description  of  the  governing  equations  or  algorithms  used  in  CCFM. VENTS  and 
their  technical  basis.  This  was  done  in  Sections  2,  where  the  generic  CCFM 
equation  set  was  presented  and  in  Section  3 where  the  detailed  technical  basis 
of  the  CCFM.VENTS-specific  equations  were  presented.  Section  4 presented  an 
overview  of  the  task  of  specifying  a particular  CCFM. VENTS  simulation. 

Methods  of  using  CCFM. VENTS  to  simulate  the  phenomena  of  wind  loading  and 
stack  effect  were  also  discussed  in  this  section. 

The  other  three  parts  of  this  CCFM. VENTS  documentation  are:  Part  II:  Software 
Reference  Guide  [1.1.1];  Part  III:  Catalog  of  Algorithms  and  Subroutines 
[1.1.2];  and  Part  IV:  User  Reference  Guide  [1.1.3]. 

Part  II  presents  the  generic  and  CCFM.VENTS-specific  features  of  the  framework 
of  CCFM  software.  This  includes  features  of  both  program  and  data  structures. 
Numerical  considerations  used  to  treat  the  solution  to  the  CCFM  equation  set 
are  also  presented.  Part  II  concludes  with  a presentation  of  methods  of  using 
the  CCFM  structure  to  generate  stages  of  the  computer  model  beyond  CCFM. VENTS. 

Part  III  is  a catalog  of  all  the  modular  algorithms  and  associated  computer 
subroutines  used  to  simulate  the  physical  phenomena  in  CCFM, VENTS.  The 
catalog  has  two  kinds  of  entries,  algorithms  which  involve  physical  phenomena 
and  other  "utility"  algorithms  required  in  the  physical  algorithm 
calculations.  Each  physical  algorithm  entry  of  the  catalog  includes  a brief 
description  of  the  phenomenon  simulated,  a concise  but  complete  presentation 
of  the  calculation  procedure  used,  identification  of  all  input  and  output 
parameters,  and  a listing  of  the  computer  coding.  These  catalog  entries 
represent  the  bridge  between  the  governing  equations  of  CCFM. VENTS  and  their 
technical  basis,  as  presented  in  Section  3 of  this  report,  and  the 
implementation  of  these  mathematical  modeling  results  into  the  software  of 
CCFM,  or  of  any  other  modular,  zone -type,  compartment  fire  model  computer 
code . 

Part  IV  is  the  CCFM. VENTS  user's  guide. 
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NOMENCLATURE 


A 

a constant,  Eq . (3.3.24) 

B 

a constant,  Eq . (3.3.24) 

b 

width  of  a vent 

C 

vent  flow  coefficient,  Eq . (3.3.14) 

C 

C at  y 

Cd 

flow  coefficient  for  entering  duct  flow 

*^INCOMPRESS  IBLE 

C for  incompressible  flow 

^PRESSURE 

local  pressure  coefficient 

Cp 

specific  heat  at  constant  pressure 

Cv 

specific  heat  at  constant  volume 

, PEN 

concentration  of  product  of  combustion  k in  a flow 
slab  as  it  penetrates  the  receiving  room 

“^k  , U [ ^k  , L ] 

concentration  of  product  of  combustion  k per  unit  mass 
of  upper  [lower]  layer  of  a room 

‘^k  , U , i t ^k  , L , i ^ 

concentration  of  product  of  combustion  k per  unit  mass 
of  upper  [lower]  layer  of  room  i 

^k  , 1 t ^k  , 2 ^ 

concentration  of  product  of  combustion  k in  room  1 
[room  2]  at  the  elevation  of  a vent  flow  slab 

^0  2 , U [ ^0  2 , L ] 

concentration  of  O2  per  unit  mass  of  upper  [lower] 
layer  of  a room 

*^02  , U , i [^02  , L , i ^ 

concentration  of  O2  per  unit  mass  of  upper  [lower] 
layer  of  room  i 

^02,1  ^0  2 , 2 ] 

concentration  of  O2  in  room  1 [2]  at  the  elevation  of 
a vent  flow  slab 

^02  , PEN 

concentration  of  O2  in  a flow  slab  as  it  penetrates 
the  receiving  room 

F 

a constant,  Eq . (3.2.3) 

Fu  [FJ 

fraction  of  flow  slab  deposited  into  upper  [lower] 
layer  of  receiving  room 
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f 

vent  flow  function,  Eq . (3.3.13),  or  fan  curve 
function,  Eq . (3.4.23) 

factorj^ 

multiplying  factors  to  predict  Pj^  source  Q,  Eq . 

(4.3.1) 

g 

acceleration  of  gravity 

Hi 

depth  of  vents  at  the  endpoints  of  a fan/duct  system 
in  room  i 

Hu,i  [Hl_J 

magnitude  of  the  difference  in  elevation  in  room  i 
between  the  top  [bottom]  of  the  vent  and  the  layer 
interface 

I R 

indicator  of  flow  direction  of  a vent  flow  slab 

JTERM 

number  of  points  used  to  approximate  the  fan  curve  > 

Hu.i 

mass  in  upper  [lower]  layer  of  room  i 

M 

^^PLUME 

mass  flow  rate  in  the  plume  as  it  penetrates  the  layer 
interface 

M 

“duct 

mass  flow  rate  through  the  duct 

Hduct.i  [HuucT.a^ 

mass  flow  rate  through  fan/duct  system  if  its  source 
is  room  1 [room  2] 

M 

“slab 

mass  flow  rate  in  a flow  slab 

HsLABU,iFR  [HgLABL, 

^Pj^]  mass  flow  rate  in  a flow  slab  entering  a fan/duct 
system  from  the  upper  [lower]  layer  of  room  iFR 

Hu.i  [ML.i] 

sum  of  mass  flow  rate  contributions  to  the  upper 
[lower]  layer  of  room  i from  all  plumes,  jets, 
sources,  etc. 

Mu  [MJ 

a mass  flow  rate  to  the  upper  [lower]  layer 

N 

^^PROD 

number  of  products  of  combustion  in  the  simulation 

N 

^^R  0 0 M 

number  of  rooms  in  the  simulation 

N 

‘•^SLAB 

number  of  flow  slabs  through  a vent 

N 

V E L E V 

number  of  unique  elevations  of  the  endpoints  and  zeros 
of  linear  segments  of  a Ap^  ^ distribution  along 
entire  height  of  a vent 
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N 

E L E V 

number  of  unique  elevations  of  the  endpoints  and  zeros 
of  linear  segments  of  a Ap^  ^ distribution  along 
entire  length  of  a common  wall  segment 

, U , i [ , 

L amount  of  product  of  combustion  k in  upper  [lower] 

layer  of  room  i 

^k  , S L A B 

flow  rate  of  product  of  combustion  k in  a flow  slab 

^k  , SOURCE 

fire  generation  rate  of  product  of-  combustion  k 

^k  , U , i [ ^k  , 

L ^ ] sum  of  product  of  combustion  k flow  rate  contributions 

to  the  upper  [lower]  layer  of  room  i from  all  plumes, 
jets,  sources,  etc. 

^k  , U [ ^k  , L 

] a product  k flow  rate  to  the  upper  [lower]  layer 

^02 , SLAB 

O2  flow  rate  in  a flow  slab 

p 

■^02  , SOURCE 

rate  of  O2  generated  by  the  fire 

P 

absolute  pressure 

Pd  ATUM 

absolute  datum  pressure 

Pn 

P at  y^, 

Pi 

absolute  pressure  in  room  i 

PfLOOR , i 

p at  elevation  yFL00R,i  in  room  i 

Q* 

dimensionless  Qsource > • (3.2.7) 

Q 

total  energy- release  rate  of  the  fire 

Qslab 

enthalpy  flow  rate  in  a flow  slab 

^SOURCE 

(1  - 

QtOT AL 

(1  - At)Q 

Qu  , i [ Ql  , i 

] sum  of  enthalpy  plus  heat  transfer  plus  energy  release 

flow  rates  to  the  upper  [lower]  layer  of  room  i from 
all  plumes,  jets,  sources  etc. 

Qu  [Ql] 

a flow  rate  of  enthalpy,  heat  transfer,  and/or  energy 
release  to  the  upper  [lower]  layer  of  a room 

R 

gas  constant 

Ri  [R2] 

duct  flow  resistance  for  forward  [backward]  flow 
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S. 

indicator  of  flow  direction,  Eq . (3.4.2) 

T 

•^D  UC  T 

absolute  temperature  of  flow  in  the  duct 

T 

•^PEN 

absolute  temperature  of  a flow  slab  as  it  penetrates 
the  receiving  room 

T 

•^PLUME 

average  temperature  of  fire  plume  gases  at  layer 
interface 

Tr 

absolute  temperature  of  the  layer  in  the  receiving 
room  which  includes  the  elevation  ygLAB 

Tu  [TJ 

absolute  temperature  of  upper  [lower]  layer  in  a room 

Tu,x  [TL,i] 

absolute  temperature  of  upper  [lower]  layer  in  room  i 

Ti  [T^] 

absolute  temperature  in  room  1 [room  2]  at  the 
elevation  of  a vent  flow  slab 

Vi 

volume  of  the  room  i 

Vu,i 

volume  of  upper  layer  in  room  i 

Vduct 

volume  flow  rate  through  a fan/duct  system 

VdUCT,1  fV^uCT, 

2 ] volume  flow  rate  through  fan/duct  system  if  its  source 

is  room  1 [room  2] 

Vp 

volume  flow  rate  through  the  fan 

Vp,  j 

volume  flow  rate  corresponding  to  point  J on  the  fan 
curve 

V 

velocity  in  a flow  slab  as  it  enters  receiving  room 

''^W  I ND 

wind  speed 

W 

Eq.  (3.3.12) 

W 

w at  y 

X 

1 - e 

X 

X at  y 

y 

arbitrary  elevation  within  a room  of  interest  above 
datum  elevation 

y 

is  an  arbitrary  fixed  value  of  y within  a room 

y 

mid-elevation  of  flow  slab 
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yCE  I L 

elevation  of  ceiling  of  a room  above  datum  elevation 

yc  E I L , i 

elevation  of  ceiling  of  room  i above  datum  elevation 

y?  I RE 

elevation  of  fire  above  datum  elevation 

yp  LOOK 

elevation  of  floor  of  a room  above  datum  elevation 

yp  LOOK , i 

elevation  of  floor  of  room  i above  datum  elevation 

yL  A YER 

elevation  of  layer  interface  of  a room  above  datum 
elevation 

yLAYER  , i 

elevation  of  layer  interface  of  room  i above  datum 
elevation 

yMIN  [yMAX ] 

minimum  [maximum]  bounds  for  y of  a vertical  vent 
between  two  adjacent  rooms,  Eq.  (3.3.5) 

yN 

elevation  N,  above  the  datum  elevation,  where  a Ap^  ^ 
distribution  has  an  endpoint  or  a zero 

yREP 

reference  elevation  above  datum  elevation 

ySL  AB 

characteristic  elevation  of  a flow  slab 

ySYSTEM  , i 

elevation  above  datum  elevation  of  the  midpoints  of 
the  inlet  and  outlet  end-points  of  the  fan/duct  system 
in  room  i 

yVTOP  [yVBOT ] 

elevation,  above  a datum  elevation,  of  top  [bottom]  of 
a vent 

z 

elevation  of  layer  interface  above  fire 

yv  T 0 P [ yv  B 0 T ] 

elevation,  above  a datum  elevation,  of  top  [bottom]  of 
a vent 

7 

C /C 

p / V 

Ap 

^Pl.2 

'^Pad  J , 1 I^^Pad  j , 

2]  adjusted  pressure  difference  for  fan/duct  flow  from 

room  1 [room  2],  Eq . (3.4.25) 

^Pduct 

pressure  drop  across  the  duct,  Eq . (3.4.2) 

^Pentry 

pressure  drop  across  the  duct  entrance  region 
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App 

cross-fan  pressure  difference,  Eq . (3.4.1) 

^Pf  , j 

cross-fan  pressure  difference  corresponding  to  point  J 
on  the  fan  curve 

^Pf , soln 

value  of  ApF  for  a valid  solution  to  the  fan/duct  flow 
problem 

^Pf  , 0 

solution  to  f(App  q)  = 0,  Eq . (3.4.24) 

^Pn 

Ap  at  Yj, 

^PSYSTEM 

'^PsYSTEM.l  ' ^PsYSTEM,2 

Api  2 

Pj^  - P2  at  a specified  y 

^^EPS 

minimum  temperature  difference  required  for  layer 
initiation 

AT 

^■^PLUME 

1 "^P  LOME  ‘ 1 

Ap 

Pi  - Pz 

6 

Eq.  (2.2.1) 

^Pd  u c t 

pressure,  relative  to  the  datum  pressure  at  a location 
somewhat  inside  the  entry  region  of  a forced 
ventilation  duct 

^PfLOOR , i 

Pfloor  , i ■ Pdatum 

'^Pp  , 1 [ '^Pp  , 2 ] 

pressure,  relative  to  the  datum  pressure,  measured  at 
fan  inlet  [outlet]  section 

5p, 

Pi  * Pfloor , i 

^Pquiescent  [^PwIND^ 

outside  pressure,  relative  to  the  datum  pressure, 
local  to  the  vent  in  a quiescent  [wind]  environment 

^Pref 

pressure  of  ambient,  relative  to  datum  pressure,  at 
elevation  y^gp 

'^Prefnew 

new  value  for  5pj^pp  , Eq . (4.5.2) 

^PSYSTEM,!  t^PsYSTEM, 

2]  pressure,  relative  to  the  datum  pressure,  measured  at 
fan/duct  system  inlet  [outlet]  endpoint 

€ 

Ap/Pi 

e 

€ at  y 

A 

a constant,  Eq . (3.3.23) 
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effective  fraction  of  Q radiated  from  the.  combustion 
zone  and  the  plume 

total  fraction  of  Q lost  to  the  bounding  surfaces  of 
fire  room 

fraction  of  a penetrating  flow  slab's  "enthalpy  of 
buoyancy"  [see  Eq . (3.3.42)]  transferred  to  the 
surfaces  of  the  receiving  room 

p 

density  of  the  layer  at  y 

^AMB 

density  of  ambient 

^DUCT 

density  of  flow  in  the  duct 

^DUCT , 1 ^^DUCT 

2 ] density  of  flow  in  the  duct  if  its  source  is  room  1 

[room  2] 

^PEN 

density  of  a flow  slab  as  it  penetrates  the  receiving 
room 

Pm  [Pr] 

density  of  upper  [lower]  layer  in  a room 

PM.i  [pL,i] 

density  of  upper  [lower]  layer  in  room  i 

Pi  [P2] 

density  in  room  1 [room  2]  at  the  elevation  of  a vent 
flow  slab 

subscripts 


iREC 

index  of  room  receiving  a flow  slab 

iFR 

index  of  room  from  which  fan/duct  system  flow  is 
extracted 

iTO 

index  of  room  to  which  fan/duct  system  flow  is 
delivered 
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APPENDIX  - Derivation  of  the  Generic  Equation  Set  for  An  Arbitrary  Room 

Consider  an  arbitrary  room  i of  a simulated  facility.  Define  as  the  fixed 
volume  of  the  room  and  Vy  ^ as  the  changing  volume  of  its  upper  layer,  i.e. , 
i volume  of  the  lower  layer.  Also,  consistent  with  the 

definitions  in  Section  2.1,  My  My  ^ and  Pk  u i > l i defined  as  the 

total  Instantaneous  mass  and  amount  of  product  of  combustion  k in  the  upper 
and  lower  layers,  respectively. 

In  this  Appendix,  the  generic  equation  set  for  the  room  will  be  derived  in 
terms  of  the  solution  variables  My  ^ , My  ^ u i > l i ’ i ’ 

^Pfloor , i ■ 


Conservation  of  Mass: 


Conservation  of  mass  in  the  upper  and  lower  layers  requires  that 


dMy^ 

y/dt 

= My,y 

(A. 

1) 

dMy^ 

y/dt 

= 

(A. 

CM 

Consistent  with  the  definitions  in  Section  2.1,  the  terms  My  ^ and  My  ^ on  the 
right  hand  side  of  Eqs . (A.l)  and  (A. 2)  each  represent  the  sum  of  individual 
mass  flow  rate  contributions  to  the  layers  from  all  plumes,  jets,  sources, 
etc.  (positive  values  mean  net  rate  of  flow  is  into  the  layer).  The  terms 
^k  u i ’ ^k  L i ’ Qu  i > i > to  be  introduced  below,  will  represent 

analogous  product  and  enthalpy  flow  rate  contributions,  respectively. 


Conservation  of  Energy: 


Conservation  of  energy  in  the  upper  and  lower  layers  requires  that 

, i , i ^ + PfLOOR  , i , i ^ Qu  , i (A.  3) 

d [My  y C^Ty  i ] /dt  * P f L 0 0 R , 1 , 1 Z ^ ^ , 1 ( A- . 4 ) 

where  C^  is  the  specific  heat  at  constant  volume,  assumed  here  to  always  have 
a constant  value  associated  with  a perfect-gas  model  of  air;  Pfloor  i 
assumed-uniform  absolute  hydrostatic  pressure  at  the  floor  elevation  of  the 
room;  and,  with  the  use  of  the  approximation  of  Eq.  (2.2.3),  the  right-hand 
terms  on  the  left-hand  sides  of  Eqs.  (A. 3)  and  (A. 4)  represent  the  rate  at 
which  pressure  forces  do  work  on  the  moving  layer  interface. 
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Conservation  of  Products  of  Combustion : 


Conservation  of  product  of  combustion  k in  the  upper  and  lower  layers  requires 


. u , i (A..  5) 

dPk.L.i/dt  = Pk.L.i  (A. 6) 


Equation  of  State  and  the  Laver  Densities.  Temperatures,  and  Concentrations  of 

Products  of  Combustion : 


Applying  uniformly  the  equation  of  state  for  air  and  using  again  Eq . (2.2.3) 
leads  to 


PfLOOR,!  , i , i ^ 


(A. 7) 


where  R is  the  gas  constant  for  air  and  Py  ^ and  p^  ^ , the  density  of  the 
upper  and  lower  layer,  respectively  would  be  computed  according  to  their 
definitions  from  the  solution  variables  My  ^ , My  ^ , and  Vy  ^ (an  explicit 
equation  for  Vy  ^ is  presented  below) . 


Pu,  i 

(A. 8) 

Pl,  i 

= - ^U.i) 

(A. 9) 

Having  computed  the  densities,  the  absolute  temperatures  of  the  upper  and 
lower  layer,  Ty  ^ and  Ty  ^ , respectively,  can  be  computed  from  the  solution 
variable  Pfloor  'i  '^sing  Eq . (A.  7). 


'^U  , i PfLOOR  , i/^Pu  , (A. 10) 

, i ^ PfLOOR , i /(Pl , i (A. 11) 


Finally,  Cj^  y ^ and  y . , the  concentration  of  product  k per  unit  mass  of 
the  upper  and  lower  layer,  respectively,  is  related  to  the  solution  variables 


Mu.i.  My,, 

P 

’ ^k  , U , i > 

^k  L i according  to 

^k  , U , 

i ""  , U , 

(A. 12) 

, L , 

i ""  ^k  , L , 

i/Ml.x 

(A. 13) 
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Deriving  Equations  for  the  Solution  Variables  Pp l q o R i 


•U  , i AYER,  i • 


and  V 


Using  Eq . (A. 7)  with  Eqs . (A. 8)  and  (A. 9)  to  replace  ^ ^ and  ^ ^ in 

Eqs . (A. 3)  and  (A. 4)  leads  to 

(Cv/R)  h [Pf  looR  , i , i ^ + Pp  L OOR  , i *^^0  , i , i (A. 3') 

(C^/R)d[ppLQOR  , i (^i  ■ ■ PpLOOR  , i , i ^ *^L  , i (A. 4') 

Now  add  Eqs.  (A. 3')  and  (A. 4')  and  use  the  result  to  replace  Eq . (A. 4') 


dPpLOOR.i/dt  = [R/(C,V,)](Qu 


Qt.i) 


Expanding  Eq . (A. 3')  leads  to 


(A. 4' ' ) 


(C,^/R)Vy  ^ ^ dppLooR  , i ■'■(^v/^^PfLOOR  , i , i PpLOOR  , i , i , i 


Collecting  coefficients  of  dVy  ^/dt  in  the  above  and  using  Eq . (A.4'')  results 
in  the  following  equation  for  ^ which  is  used  to  replace  Eq . (A. 3'): 

dV,_.dt  = [(1  - - (Vu,iZV,)QL,  J/[PpL00R.i(l  + c^ZR)] 

(A. 3") 


Introducing  7,  the  ratio  of  specific  heats, 

7 - Cp/C,  = 1 + R/C^  (A. 14) 

Eqs.  (A.3'')  and  (A.4'')  are  written  finally  as 

dV^  ,/dt  = [(7-l)/(7PpL00R,i)][(l  - - (Vu.iZVi)QL.i] 

(A.3"  ' ) 

dPpLOOR.i/dt  = [(7-1)ZVJ(Qu,,  + QL.i)  (A.4'") 
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For  rooms  i of  a facility  with  horizontal  sectional  areas,  , which  are 
uniform  with  elevation,  (i.e.,  for  rectangular  rooms  with  vertical  walls  and 
for  other  right  cylindrical  rooms),  Eq . (A.3''')  can  be  written  in  terms  the 
solution  variable  Ylayer  i ’ elevation  of  the  interface  separating  the 

upper  and  lower  layers,  instead  of  Thus,  if 


(yCE  I L I NG  , i ‘ yFLOOR.i)  (A. 15) 

, i ^ E I L I NG  , i ’ yLAYER.i)  (A. 16) 

where  Yceiling  i ypLOOR  i fixed  elevations  of  the  ceiling  and 

floor,  respectivelY , then  Eq . (A.3''')  can  be  replaced  bY 

^yt  AYER  , i ■ ^^/^yPpLOOR  , ^ H (ycE  I L I NG  . i ’ yL  A Y E R . i ^ , i 

^yLAYER,i  ' yp L OOR , i , i ^ 

(A.3"  " ) 


The  Generic  Equation  Set: 


The  final  set  of  generic  equations  for  the  solution  variables  My  ^ , M^  ^ , 

^U.i-  Pfloor.i-  Pk.u.i-  is  Eqs.  (A.l),  (A. 2),  (A.3'"),  (A. 4'"), 

(A. 5),  and  (A. 6),  respectivelY-  Also,  the  derived  variables  py 

Pi  i>  i’  '^L  i>  u i’  ‘^k  L i Computed  from  Eqs.  (A . 8 ) - (A . 13 ) . 


For  right  cYlindrical  rooms  where  Eqs.  (A. 15) 
Variable  Ylayer  i used  instead  of  the 
(A.3'  ' ')  is  replaced  bY  Eq . (A.3"  ")  . 


and  (A. 16) 
variable  Vy 


are  satisfied,  the 
^ , in  which  case  Eq . 


The  Perturbation- Pressure  Variable : 


TYpical  fire  scenarios  involve  PpLooR  i which  varY  a relativelY  small 

amount  from  a datum  absolute  pressure,  Pdatum  • example,  if  Pdatum  ^ 

problem  is  taken  as  the  local  atmospheric  pressure  at  ground  elevation,  of  the 
order  of  10^  pascals,  then  tYpical  variations  of  the  Pplqor  ftom  Pdatum 
will  onlY  be  of  the  order  of  10^  or  less.  For  this  reason,  a new  transformed 
pressure  variable. 


5p 


floor , i 


Pf 


LOOR , i 


- Pd 


ATUM 


(A. 17) 
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will  lead  typically  to  a higher  level  of  accuracy  in  the  computation  of 
Pfloor  i • perturbation-pressure  variable  SppLooR  i adopted  as  a new 

solution  variable  to  replace  Pflqor  i • 

The  generic  equation  set  using  the  perturbation-pressure  is  presented  in 
Section  2.3. 
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The  two -layer  environment  for  an  arbitrary  room  of  a multi -room 


Figure  2.1. 
facility . 
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FLOOR,! 

vxvs^gtum  elevation 


Figure  3.2.1.  A room  of  fire  origin. 


ICEIL,! 


wXXX 
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»))))))))))))))))))))) 

Configuration  4 

(d) 


(b) 


(e) 


Configuration  6 
(f) 


Figure  3.2.2.  Configurations  for  the  fire,  the  plume  and  the  elevations  of 
interface  and  ceiling. 
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Figure  3.3.1. 


The  fire -generated  environment  in  two  adjacent  rooms. 
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Figure  3.3.2.  Possible  distributions  of  2 ^ Pi  ‘ P2  ^ common-wall 

segment. 
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Figure  3.3.3.  Possible  distributions  of  Ap^  2 ^ Pi  ' P2  ^ single  vent  and 

the  flow  directions  of  its  uniform-property  flow  slabs. 
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Conditions  in  adjacent 
pair  of  rooms 


Pressures  and  pressure 
drops  across  the  common 
wall  segment 


Figure  3.3.4.  Use  of  the  algorithms  DELP,  COMWLl , VENTHP,  and  PX0G02  to  carry 
out  rates  of  exchange  of  mass,  enthalpy,  oxygen,  and  other  products  of 
combustion  through  vertical  natural  vents  of  a multi -room  facility. 
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Figure  3.4.1.  The  fire -generated  environment  in  two  rooms  connected  by  a 
fan/duct  system. 
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Figure  3.4.2.  A sketch  of  a generic  fan  operating  curve  which  increases 
monotonically  in  volume  flow  rate  with  decreasing  cross -fan  pressure 
difference . 
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Figure  3.4.3.  An  illustration  of  the  flow  entering  the  duct  in  room  IFR  when 
the  layer  interface  is  between  the  top  and  bottom  of  the  vent  opening. 
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Ap 


FAN  CURVE  : 


Figure  3.4.4.  Graphical  presentation  of  the  problem  of  solving  Eq.  (3.4.24) 
and  finding  the  direction  and  flow  rate  through  a fan/duct  system. 
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Figure  3.4.5.  Graphical  presentation  from  Figure  3.4.4  of  the  four  classes  of 
solution  for  flow  rate  through  the  fan/duct  system:  (a)  single  forward-flow 
solution  and  no  backward- flow  solution,  (b)  single  backward- flow  solution  and 
no  forward-flow  solution,  (c)  one  forward-flow  solution  and  one  backward-flow 
solution,  and  (d)  no  solutions. 
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Conditions  at  two 
end-points  of  fan/duct 
system  and 

characteristics  of  system 


VENTF 


Direction  and  properties 
of  fan/duct  flow  as  it 
enters  and  exits  fan/duct 
system 


If  flow  is 
unspecified 


Forced  vent 

contributions  to  rates  of 
mass,  enthalpy,  and 
products  flowing  to  the 
upper  and  lower  layers 
of  rooms  joined  by  the 
fan/duct  system 

I 


Figure  3.4.6.  Use  of  the  algorithms  VENTF,  FANRES , and  FL0G02  to  carry  out 
rates  of  exchange  of  mass,  enthalpy,  oxygen,  and  other  products  of  combustion 
through  fan/duct  forced  ventilation  systems  connecting  two  rooms  of  a multi- 
room facility. 
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